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A STUDY OF ZODIACAL LIGHT MODELS 


SUMMARY 

Various interplanetary dust cloud models which were formulated to 
predict the photometric properties of the zodiacal light are considered. A 
compilation and analysis of the models in the literature is presented. This 
study revealed a lack of unified coverage of the zodiacal light models in the 
literature . 

The important points of the present study are as follows: 

1 . A detailed compilation of most of the models presented in the 
literature. 

2. A correlation of the models and their association with the develop- 
ment of more detail models . 

3. The coverage of out -of -the -ecliptic models. 

4. An introduction to the inversion formulation. 

The authors' own calculations of single particle Mie scattering results, size 
distribution Mie scattering results, and ecliptic and out-of-the-ecliptic zodi- 
acal light model results are presented to illustrate and extend the presentation 
for various models. 

The conclusion of the study is that at present there is not a unique 
model for the zodiacal light, although extensive work has been carried out 
by various authors. However, it is believed that this direction of research 
is fruitful and, with new observation results, defensive data on the nature of 
the interplanetary dust cloud can be obtained. 

The study makes the following recommendations for better refinement 
for model calculations and analysis of the zodiacal light by photometry: 

1. The size and spatial distribution needs to be based on experimental 
data. This is particularly true in the case of particle size distributions where 
a lower cut-off limit is assumed. 



2. A set of polynomials needs to be created to form do/d£2 . 

3. The effects of irregularly shaped particles compared to spheres 
should be better defined to categorize differences. 

4. The method of inversion of the photometric observation to find 
particle sizes and composition. needs more definition and study. 

5. The complete effects of the physical forces on the dynamics of 
the interplanetary particles need to be determined for input into models. 

6. A Mie theory model which satisfies all the gross properties of 
the observation data should be found, if possible. 
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SECTION I . INTRODUCTION 


Zodiacal light is produced by sunlight scattered. by : the interplanetary 
cloud of particles around the sun. The cloud is concentrated toward the plane 
of the ecliptic so that the zodiacal light appears as a band along the ecliptic. 

In the vicinity of the sun, the interplanetary particles merge with the .plasma 
surrounding the sun to form part of the solar corona. In the antisolar direction, 
a slight increase in brightness occurs; this phenomenon of the zodiacal light is 
called either gegenschein or counterglow. ’« * * : 

Zodiacal ‘light' observations have a lorig history. . It;is known that zodi- 
acal light was observed by the Egyptians and- the Greeks [106] . The. first 
recorded systematic series of observations of zodiacal light was made by - 
Cassini in 1672 [20] . Cassini proposed that the zodiacal light was sunlight 
scattered by. an interplanetary lenticular dust cloud. Searle [ l] compiled the 
notes on the observations of the zodiacal light and gegenschein made at the 
Harvard University Astronomical Observatory during the period from 1840 to 
1890. He modeled the zodiacal cloud using different phase functions and spa- 
tial densities. He found that the best fit to observational data was given by a 
phase function approximately proportional to the phase angle and a density 
inversely proportional to the distance of the scatterers from the sun, although 
data from asteroids indicated to him that it may fall off faster. To account 
for the increase in brightness near the sun, an increase in the amount of 
meteoric matter was suggested. Fessenkoff in 1914 attempted further theo- 
retical models to explain his own observations [106]. He assumed the meteor 
cloud to be infinite in extent and the density to vary inversely with the distance 
from the sun. For phase functions, he used Euler, Lambert, and Lommel- • 
Seeliger laws. 

As reported by Schoenberg [2] , Van Rhijn [3] , in the course of at- 
tempting to determine the contribution of scattered stellar light to the sky 
brightness, observed that zodiacal light constituted almost 43 percent of the 
light in the night sky. Van Rhijn modeled the zodiacal cloud, assuming phase 
functions for the particles that corresponded to meteor phase functions and 
assuming that the spatial density was constant over spheroids. His studies 
and calculations extended outside the ecliptic. Dufay [4] reported a substan- 
tial polarization of the zodiacal light through photographic analysis. Elvey 
and Roach [5] made detailed photoelectric observations. Whipple and 
Grossner in 1949 [ 107] used the polarization data of Dufay and the brightness 
data of Hoffmeister [6,7] and Elvey and Roach [5] to set upper limits on the 
electron density. However, they concluded that the dust component alone could 
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probably account for the observed polarization. Behr and Siedentopf [8] made 
the first accurate measurements of the polarization. From 'these measurements 
it was proposed that near the sun (e < 35 deg) , almost one-half of the zodiacal 
light was because of scattering from electrons. _ , i 

The spectrum of the zodiacal light was observed by Fath [108] in 1909 
and found to be like that of the sun. Blackwell and Ingham [9, 10, 11] used the 
spectrum of the zodiacal light obtained by them to infer an upper limit of the 
electron density much lower than that" predicted by Behr and Seidentopf. 

Alien [12] and Van de Hulst [13] independently showed that the increase 
in brightness of the zodiacal light near the sun was caused by light diffracted by 
the interplanetary particles. Because of this result the requirement for an 
increase in spatial density of the particles near the sun could be omitted. 

Richter [ 14] made studies of the scattering of light from dust clouds 
using scattering functions measured in the laboratory. Walter [15] and Giese 
and Siedentopf [16] used the Mie theory of scattering from spheres to obtain 
theoretical scattering functions for use in gegenschein and zodiacal light model 
calculations . 

Detailed radiance and polarization measurements of the zodiacal light 
were made by Weinberg [17, 18] . Smith, et al. [19] obtained a whole sky 
radiance map. 

It is impossible to include all the research in a condensed history of 
the zodiacal light as is attempted above and summarized in Table 1. Such 
abbreviation inevitably omits the works of many important contributors. The 
above history is presented to give an indication of the work done in the past 
on zodiacal light. Since about 1960, there has been a great deal of study 
(theoretical and observational) of the zodiacal light, using modern tools and 
techniques, which does not show in the condensed history. 

Representative results* of photometric observations of the radiance and 
polarization in the ecliptic are shown in Figure 1. Figure la presents the 
range of the results for the radiance of the zodiacal light. The ordinate is 
radiance in units of the number of 10th magnitude stars per square degree and 
of the average brightness of the solar disk (B^) . The abscissa is the elonga- 
tion in degrees; elongation is the angle between the earth-sun line and the line 
of sight. The radiance increases rapidly as the zodiacal light merges with the 
solar corona. Gegenschein is indicated by a much smaller peak in radiance in 
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3 1200 


ELONGATION (DEGREES) 


ELONGATION (DEGREES) 


Figure 1. Representative results of (a) radiance and (b) polarization of the 
zodiacal light in the plane of the ecliptic. (The solid line indicates the average 
value, and the shaded region indicates the spread in the observational data. ) 






r TABLE 1. A CONDENSED HISTORY OF THE ZODIACAL LIGHT. 


Egyptians [106] Observed the zodiacal light and drew pictures of it .... 

; - ■ in the pyramid. «-■ : ' 

Cassini [20] , Observed the zodiacal light and proposed a lenticular 

cloud, of dust in the ecliptic for its origin. " 

■ ' ;• -- . i ■ ’ " . ( 

Brorsen [109] Observed the gegenschein. 

Searle [1] Calculated the radiance by using asteroidal observations. 

Compiled observations at HCO. : ( . . . . 

Fath [108] Photographed the spectrum of the zodiacal light. 

Van Rhijn [3] Observed the zodiacal light photometrically and calculated 

dut-of-the-ecliptic models for ;the radiance of the zodiacal 
light. 

Dufay [4] Found considerable polarization in zodiacal light. 

Elvey/Roach [5] Carried out detail photoelectric observations. 

Allen [1 ! 2] arid ' Added diffractionpeak contributionto the F-corona calcu- 

Van de Hulst [13] « lations; . t . - . e : . . 

Richter; [14] Made ( laboratory studies of the scattering of dust clouds. 

Walter [1 5] Used the Mie theory in gegenschein model calculations. 

Blackwell/lngham Used the spectrum of the zodiacal light to set an upper 

[9,10,11] limit of 120 electrons/cm 3 at the earth in interplanetary 

■ - ■■ •■■■ spacer . , .. } 

Giese/Siedentopf [16] Used the Mie theory to obtain optical properties of th'e 

zodiacal light for various models. 

Weinberg [17, 18] Obtained detailed radiance and polarization data. 


Smith et al. [19] Obtained a whole sky radiance map of the zodiacal light. 




the antisolar direction. In Figure lb , the range of the degree of polarization 
of the zodiacal light versus elongation is presented. The maximum degree of 
polarization occurs at approximately 70-degree elongation. 

The particles which produce the zodiacal light are assumed to range 
between 0. 01 and 10 microns in diameter with a peak frequency distribution 
near 0. 5 micron. The spatial concentration is toward the ecliptic and decreases 
as a function of radial distance from the sun. This distribution is supported by 
photometric measurements and by analysis using scattering theory. 

A number of hypotheses have been generated to explain the existence Of 
interplanetary dust. These hypotheses include the following [21] ; 

1. Dust injection from comets. 

• ■ r • i • i , , . . . ■ • . ■ 

2. Debris from asteroidal collision. 

3. Grains from condensation of gas on nucieation centers. 

! ’ v • .1- : > . ... 

4. The remains of an original and very dense dust cloud. 

j 5. The passing of the solar system through an interstellar dust 
cloud. 

i 

In recent years, there has been renewal of interest in problems con- 
cerned with the interplanetary; medium, since it is believed that the determina- 
tion of the correct hypothesis for the existence of the interplanetary medium 
will be a major step forward in understanding the environment of interplanetary 
space and, by apocatastasis, the origin of the solar system. Earlier studies of 
zodiacal light were performed visually or with photographs at unfavorable sites 
and include uncertain corrections for extinction and extraneous sources. Photo- 
electric techniques are now being employed to measure the distribution over the 
celestial sphere of the parameters of the light. Earth-based observations are 
still hampered by the earth’s atmosphere. However, polarimetric observations 
from spacecraft are now planned. These observations will give better data on 
the radiance, degree of polarization, color, and orientation of the polarization 
plane of zodiacal light. 

To deduce physical characteristics of zodiacal dust from photometric 
observations, 1 it is necessary to compare the observed properties with those 
calculated for several models of the zodiacal cloud. These several models - 
must be used, for there is no unique model for the cloud because of the number 
of parameters in each model and the lack of sufficiently accurate experimental 
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data; nonetheless, the models are very informative. The difficulty in obtain- 
ing a unique model may be appreciated by considering the following major steps 
in the calculation of a zodiacal light cloud model [22] : 

1. Select a particle size and shape distribution which is consistent 
with the dynamics that act upon the particles and their physical constraints. 

2. Calculate the individual scattering functions for the particles 
assumed. 

3. Calculate the photometric properties, based on the spatial distri- 
bution assumed for the model. 

4. Analyze the results of varying the parameter to determine the in- 
dividual uniqueness. 

The last two steps are presented in detail later in this report. 



SECTION II. SCATTER I NG THEORY 
A. Introduction 


To be able to. accurately model a cloud of scattering particles, the 
scattering properties of the individual particles composing the cloud must be 
known. In this section, the general scattering theory for individual particles 
will be presented and then extended to a tenuous cloud of such particles. 

Since the scattering is classical field scattering for electromagnetic radiation, 
the derivations are omitted since they are available in the literature. This 
section presents the general treatment for scattering, the forms of some 
differential scattering cross sections which are used to describe the angular 
distributions of the scattered field from different scatterers, and a discus- 
sion of the size distribution for the scattering centers. 

B. General Scattering Theory 

To know the scattering properties of the particles, one must know the 
relation between the electromagnetic fields that are incident on and scattered 
by the particle. Such a relation depends on the shape and physical properties 
of the particle and, in general, the relation is very difficult to describe 
analytically. However, the scattered properties can be described in a general 
formalism without direct reference to the properties of the particle. Con- 
sider the electromagnetic fields at a large distance r from the particle. The 
field is transverse to the direction of propagation, and the amplitude of the 
scattered waves, is a function of position. The fields are a combination of the 
incident wave, considered hereafter to be a plane wave with a wave vector 
k= 27 t/X, and a spherical outgoing wave. 


EVE 


o 


e i ~~ r + A (0, <p) e 1 


( 1 ) 


where I A (0, 0) | is the amplitude of the scattered wave, e^ r is its phase, 
and |-is the geometric attenuation factor. (0, 0) is such that the scattered 
part contributes nothing to E before scattering. It is related to the incident 
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field by a scattering matrix, S, which will be'defined shortly. The linear 

relationship between the incident and scattered waves is a consequence of the 
linearity of Maxwell's equations and is valid for the cases considered here. 
Before defining S, consider Figure 2 which shows the coordinate system used 

to describe both the direction and polarization of the incident and scattered fields 



Figure 2. The geometry for incident and scattered;, radiation, 
and the scattering plane from a different perspective. 

The X, Y, and Z axes are fixed to the scattering center and oriented 
such that the Z-axis is along the direction of the incident radiation. The unit 

vector n and the propagation vector k are in the direction of the incident radi 

1 1 

ation. The unit vector ft' and propagation vector k' are defined similarly 


E 




for the scattered radiation, n and n define the scattering plane and the 
scattering angle, 0, which is the angle between them (0 < 0 < 180 deg) . 

a ! ‘ a ■*' ■ • . ■ - ■ . •••■•> ■- 

A and A. are the same and are perpendicular to the scattering plane which is 

inclined an angle, cp , to, the X-axis. Since the amplitudes are transverse to > 

A A" A A A" A > 

k , the unit vectors 4 and 4 are perpendicular to A, n and a, n , respectively. 

The quantities and are the electric field components of the incident 

radiation A and 4 . Likewise, and are, the electric field components of 

the scattered radiation along A. and 4 . The amplitudes of the incident and 
scattered waves .are related by the scattering matrix S , 


■ i" i 1 



In general, scattering amplitude functions S , S , S , and S are dependent 

1 J O 4 

on the particle properties, complex, and functions of 0 and <p . The compo- 
nents (E., E ) and (E E ) are referenced to coordinate systems 

^ ^ ^ A* Ox O/l* 

(4 , a ) and(f', A), respectively. Hence, the scattering niatrix is a function 
of k and k *, . ••• * 

S (e,cp) = S (k , k ') : >'- ■ • 

— — r*j -i 

• l 

where S reduces to the unit matrix for no scattering. The S-matrix contains 
all the information one can obtain about the particle from a scattering analysis. 

!■ For. a spherically symmetric particle, the scattered ' component 

is directly related to the incident 4 component only, and similarly for the a 
component. Hence, the S 3 and S 4 elements of the S matrix are zero. 



In the case of a spherical particle, the energy at a given wavelength which 
crosses a unit area in unit time at a distance r is given by 




k" 2 r " 2 


*(«.*) 

5 o\ 



( 4 ) 


where 


is the spectral radiant density and is related to E by 



e is the electric permitivity and jx is the magnetic permeability. 


The scattering intensity functions 4 and i 2 are defined as the square 
of |S 1 | and |S 2 |, respectively, i t and i 2 are calculated for spherical par- 
ticles by the use of the Mie theory [23, 110] . As will be seen, i 4 and i 2 play 
an important role in modeling zodiacal light. 

Particles which have no three-dimensional rotational symmetry have 
elliptical polarization; i. e. ^ the incident £ and £ component will be scattered 
such that each will have a V and h. ’ component. Therefore, S 3 and S 4 do 
not vanish; hence the general form for i t and i 2 is 


it = I Si(0 , <p) + S 4 (0 , 0) | 2 
and 

i 2 = IS 2 (0,0) + 83 ( 0 , 0 ) | 2 

C. Differential Scattering Cross Section 

The angular distribution of the scattered radiation can be described by 
the differential cross section. The differential cross section dc/dfli is 
defined as the ratio of the flux scattered from an object in the direction 9 
and <p per unit solid angle to the flux per unit area incident on the object, 
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M-454 



da _ energy scattered/unit tLme/unit solid angle 
dfi energy incident/unit time/unit area 


In terms of the Poynting vector S for the incident, and scattered fluxes 


da 

dO 


l S 

~sca 


~inc 
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where r is the distance of the observation point from the scattering center.. 

The differential cross section has, the dimensions of area per steradian. In 

d cr . ; 

terms of the flux density — is 

d£2 



Hence, in terms of the components i^ and i 2 , the differential cross-section 
components are 



A, 


i 




The scattered flux density for unpolarized incident light is given by 


i 



Differential scattering cross sections of some well-known scatterers are 

listed in Table 2 for reference. The notation is that of Van de Hulst [120] 

except that C is polarizability. 

P* 
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■ Scatterer 

■ dcr. .... : 


dfi ..... , ... 

- l ■ \ '■ 

Isotropic Scatterer • ..> 

a 2 /4 . .. , .. 

Lambertian Sphere 

’ 2a 2 *' v ' • "• :,VI . 

■ rpr (sin0.,-0. cos 9) , •. ....... 

o7 r ’■ '■ ....... 

Circular Diffraction » ■ 

' 

2 1 2 

k a 2J (ka0) • I- ... : 

Rectangular Diffraction 

k 2 a 4 .. 

— — E(Jkb sinfi .eos^) E(Jkc sin0 cos0) 2 

Rayleigh 

C 2 J 2 

of ( 1 + cos 0) 

■ 87T€ . 1 

Electrons 

^ e 2 ^ i(l + cos 2 0) 

RayleighiGans Spheres 

fe) 2 i_lL (1 +cos 2 0) J ^ (u) 

: ' ' 8 ire . .... 

Rayleigh-Gans Cylinders 

— 4 

C P ik (1 + COS 2 0). J (u) J, (v) 
? ' 47TCU 1 5 







D. Size Distributions 


Electromagnetic radiation traveling through a scattering medium can < 
be scattered one or more times, depending on the distance between scattering j 
centers. In many situations (e.g., zodiacal light, fogs, haze, some clouds), 
the scattering centers are widely separated and, hence, the radiation is 
scattered only once before striking a detector. It is this special case of sin- 
gular scattering that is considered in this section. Another condition that is 
usually assumed in single scattering is that the scattering centers are ra.n- 
domly located with respect to each other. When both conditions are met, each- 
scattering center can be considered separately and the effects of all the centers 
is just the sum of the effects of each . ' 

The components of the differential cross, section from an isotropic, 
homogeneous cloud of particles is given 

= T 2 Z i; Cm, 0, ka) n(a). Aa ■■■;■- j = 1, 2 

k Aa 3 

where the sum is over the particle sizes. For an isotropic, heterogeneous 
cloud, the differential cross section is the sum of the cross sections of each 
particle, 

(dl). h? v (m f e - ka) n(a) Aa 

where each sum is taken only for a given refractive index (material) and all 
particles are assumed to follow the same size distribution law. 

. \ 

If the size distribution is continuous, the summation can be replaced by 
an integral 




where n (a) represents the size distribution. As can be seen, n (a) plays an 
important role in scattering. For this reason, several families of size dis- 
tributions are discussed in this section. 

There are many size distributions used in the study of light scattering; 
most are empirical formulae. Usually, two parameters are sufficient to de- 
scribe the distribution, one parameter being associated with a particular par- 
ticle size (mean, mode, median) and the other defining the spread or variance 
in the distribution. 

Distributions may be unimodel, bimodel, etc. However, since the 
simple bimodel and higher model distributions can be represented as a sum of 
unimodel distributions, only two parameter unimodel distributions will be con- 
sidered . 

The formulae for the distributions discussed below are categorized 
according to form. The four categories used are listed in Table 3. 


TABLE 3. SIZE DISTRIBUTION FORMULAE 
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Category I is a modified gamma distribution 

n t (a) = n Q a r exp(Hba P ) (8) 


which has been used by Deirmendjian (1966) to describe haze and’ clouds. 
n o» y* b, and P are real positive numbers and y is an integer. As used by 
Deirmendjian, no, y, b, and p are not independent of each other and are 
related to quantities in the frequency distribution, n is related to the total 
number of particles per unit volume, and b is relates? to the mode radius, 
provided y and p are fixed. If p = 1, nj-(a) reduces to the gamma distri- 
bution. 

Differentiating n^a) with respect to in a, one obtains 

a 3T [ V<*>1 = ? [ 1 '(t') P ] ■ (9) 

where a is the model radius. This formula represents the slope of the curve 

jf 

of log n (a) versus log a. Figure 3 shows the size distributions for some 
particular examples of n t (a) . Table 4 gives the parameters used for each 
distribution. The mode radius is the same in each one, and the parameters 
were chosen so that the slope of the curves is -4 at a radius of 1.6 microns. 

A particle density of 1.5 x 10~ 15 cm ~ 3 was used. 

TABLE 4. PARAMETERS FOR 
THE MODIFIED GAMMA DISTRIBUTION OF 
FIGURE 3. 


V 


Distribution 

y 

p 

n o 

b 

1 

6.0 

0.1705 

7.4247 x 10 +8 

54.1778 

2 

3.0 

0.2828 

3.7400 x 10 -7 

21.7002 

3 

1.0 

0.5372 

2.9998 x 10" 13 

7.2442 
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The family of distributions represented by 


\< a) 



( 10 ) 


represents Category II distributions. It includes the normal or Gaussian 
distribution 


Q o / 

n 2 (a) exp (- (a -a) 2 

'if7TCr \ 2 (T 2 


(ID 


where the parameters a and a are the mean and standard deviation 'defined by 


a* = f +aa a n (a) da V 

J -CO 2 ' 

a 2 = (a-a) 2 n (a) da 


This distribution has the disadvantage that it admits particles with ’’negative” 
radii, Figure 4 shows the Gaussian distribution for a mean radius of 1.0 
and standard deviations of 1.0 and 3.0. Only the portion of the curve corres- 
ponding ;to positive radii- is shown. Figure 5 shows the form of the Category II 
distribution for simple powers of exponent. ! 

y, t 

Category III distributions are those that are similar in form to the 
logarithmic distributions described by 


. m. . m 

hy: (a) Dgm a exp 

l f> ■ 1 

where ri is the normalizing factor for n-f 1 (a) . m is a shape factor, a is one 
Om . 

of the moments of the distribution (media, mean, mode, etc.), and' a is a 
measure of tiie width of the dispersion. Figure 6 shows a plot of equation (12) 
for several values of m. Note that this family of curves is skew symmetric. 



A simple relation between the mode radius, a, cr 2 , and the shape factor m is 
found bv differentiating: n (a) with respect to a and setting it equal to zero. 

( ina - i ria\ 

“ — — / • 
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RADIUS (a) 


Figure 4. Gaussian particle size distribution for positive radii 
(Mean radius = 1.0 andcr = 1.0 and 3.0.) 







RADIUS (a) 


Figure 5. Category II distributions for some examples of 
simple powers of the exponential. 






If m = -1, the logarithmic distributiombecomes the logarithmic normal 
distribution \ 



In this case, a_^ is both the median and the mean value of a; is the standard 

deviation for i na. The log normal distribution has the advantage that the 
median (which is also the geometric mean) remains invariant over allcr's 
(Fig. 7). 

Espenschied [111] and Kerker [24] make extensive use of the zeroth 
(order logarithmic distribution ( ZOLD) in their work on aerosols. For this 
distribution, m = 0 and equation (12) reduces to 



a is now the model value and remains invariant for all values of a. if the dis- 
\ 2 

tribution is narrow, i.e., if a is negligible compared to a, then a is related 
to the standard deviation o-g by 

; <j = a <j s ' ' . • 

As <j gets smaller," ZOLD . becomes symmetric as shown in Figure 8. 

» 

Category IV distributions are most popular for zodiacal light studies. 
The inverse power law distribution 

n 4 (a) = % (“) P -• (15) 

fj. 

represents this distribution. For published zodiacal light models, k has 
usually been in the range of 2.5 to 4.0. As with the other distributions, the 
distribution can be modified; e.g., Aller et al. [22], in an attempt to obtain a 
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Figure 8. The zeroth-order logarithmic distribution 



better model of the zodiacal light, allow k to .vary for. different ranges of a. 
This form of the Category IV distribution is the broken power law. Figure 9 
shows the Category IV distribution for the values of k noted above. It is evi- 
dent that the smaller particles dominate in this distribution; ■ • 

* • t 

The importance of the inverse power law comes from the fact that the 
cumulative mass distribution for meteors between 10 and 10”* 1 g (a radius 
between 1 cm and 1 micron) is relatively linear on a log-log plot [25] - This 
region covers the range from visual observations to dust-particle sensors on 
rockets and spacecraft. The number density increases with decreasing size, 
although there must be a physical limitation to this trend. . If and when a lower 
particle size cut-off is determined, the other categories will be more impor- 
tant in the study of. zodiacal light. 
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SECTION III. MODEL CALCULATIONS 

i • 


A. Introduction 


The main objective in the study of zodiacal light is to deduce the nature, 
amount, and distribution of the particles responsible for this light. This know- 
ledge is of fundamental importance in the theory of the origin and structure of 
the solar system. The zodiacal light can be interpreted on the basis of scat- 
tering of sunlight by a cloud of interplanetary particles. The properties of this 
cloud are to be determined without a priori knowledge of spatial extent, particle 
shape and size distribution, and concentration, and the variation of these param- 
eters with radial distance from the sun. The mathematical formulation for the 
problem is relatively simple, though there are difficulties related to the obser- 
vational uncertainties and the choice of models. 

The basic scheme is to make as many models as feasible and to approach 
a solution by trial and error, checking for consistency, and agreement with 
other theoretical criteria. To date, the number of parameters for the models 
exceeds the number of observed integral quantities. 

There are two important questions to be answered: (l) which forms of 
the mathematical functions permit the expression of the physical properties 
observed? (2) what are the values of the parameters of these functions which 
satisfy the observed quantities? It is by no means certain that there is a 
unique answer to each of these questions. The approach to the problem is 
through the integral equations which characterize the physical process of the 
scattering of the light. These equations for the radiance L are as follows: 

(e) = / / . . . / rj (a, 5 , m, x, y, z) • 1^ (a, 4, m, <p, 0) 

• da d£ dm dx dy dz d <j> d0 

j = 1, 2 . (16) 

The number density t] and the scattering intensity functions i^ are to be deter- 
mined. Generally, one postulates the form of these functions, evaluates the 
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integral equations, and compares the calculated values for the radiance compo- 
nents with observation. This approach manifests itself in the approximation of 

the functions 17 and i^ with their dependence on the particle size a, refractive 

index m, position x, y, z, and other intrinsic, variables £ . 

The calculation of the observable photometric radiance and polarization 
based on various models is presented in this section. The calculated values of 
these quantities are compared to the observed data in Section V. Radiometric 
quantities are used in the development. These quantities are defined in Table 5 

B. Cloud Radiance Integral .■ 

A basic assumption throughout the development is that of single scatter- 
ing.' For 1-micron particles with a density of 10“ 10 particles/m 3 , the mean 
free path is on the order of 10 11 AU. Single scattering appears to be an appro- 
priate assumption. Hence, extinction of the solar radiation and multiple scat- 
tered radiation will be neglected. 

• * - ; k . i 

Figure 10 shows the geometry for sunlight scattered by interplanetary, 
particles in a cylindrical volume element of length dA and volume dA. dA 
at a distance r from the sun and a distance A from the earth. 

The particles in the volume element dA g dA receive solar energy 
in the wavelength interval A to A+ 5A at the rate of 


6A [dA^ dA 77 (a, r) da] G 


(17) 


where 


/ 

G = the projected area normal to the incident sunlight (tt a 3 for 
spheres of radius a) (m 2 ) . 

= the spectral radiant flux density at the surface of the volume 
element (w/ m 2 - nm) . 

r]( a, r) = the number of particles per unit volume per unit increment 
of the size factor (radius) at a at a distance r from the 
sun (number/m 4 ) . * • . 

From the definition of the differential cross section, the fraction of energy 
received that is scattered toward the earth is 
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TABLE 5. RADIOMETRIC QUANTITIES 


Radiant Energy , Q, is the energy propagating in the form of 
electromagnetic waves (j). 

Radiant Density, W .= dQ/dV, is radiant energy per unit volume 
( J/m3) . 

Radiant Flux , 4> = dQ/dt, is the time rate of flow of radiant 
energy ( w) . •. 

Spectral Radiant Flux , = d4>/dA,'is the radiant flux per unit 
wavelength interval at wavelength A (w/.Nm) . 

Radiant Intensity , I = .d<f> /du> , of a source is the radiant flux 
proceeding from the source per unit solid angle in the direction 
considered (w/sr) . . • t 

Spectral Radiant Intensity , 1^ = dl/d^, is the radiant intensity per 
unit wavelength interval at wavelength A (w/ sr-Nm) . ■ < 

Radiance , L = d 2 4> /dco(dA cos$), in a direction, at a point on 
the surface of a source, is the quotient of the radiant flux d$ 
leaving and propagating in directions defined by an elementary 
cone containing the direction by the product of the projection of 
an element of the surface dA containing the point onto a plane 
perpendicular to the given direction and the solid angle cku of 
the cone (w/sr-Nm). 

Radiant Flux Density at a Surface (Irradiance) , E = d$/dA, is 
the quotient of radiant flux at that element of surface to the area 
of the element (w/m 2 ) . 




VOLUME 

ELEMEN 
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Figure 10. The geometry for light scattering from interplanetary medium. 


/<j£\ ^ 

\ dfi / ttst 


( 18 ) 


where co e is the solid angle the detector subtends at the volume element. For 
simplicity, spherical seatterers have been assumed. 

To obtain the radiance, the flux coming from the particles is divided 

by the area element dA and by the solid angle through which it propagates 

0 

from that area a; • The result is 
e 

zrk- < 19 > 

e e 


and 

dL = E 8\dA«(a, r) da — 
X X ' v ' dO 


( 20 ) 


dL = E 6X [dA dArj (a, r)da]7ra^ 
X X e 


(da\ ^e 
\dS2/ Tra 2 


In the literature, the term 7ra 2 1(0) may appear in lieu of dcr/dJ2 . The 
scattering function of a sphere I (0 ) is defined as the amount of radiation 
scattered per unit solid angle in a direction inclined at an angle © to the inci- 
dent direction divided by the radiation that falls on the geometrical cross 
section of the particle; 


I (0) 



( 21 ) 


Since radiance L is the flux <£ per unit area of the emitting surface 

A per unit solid angle co , then if a detector A at a distance A from the 
e e d 

surface sees the surface through a solid angle co^, one has the relation: 


L = $/A co = $/(A A /A 2 ) = $/A w ; (22) 

e e x e d ' d d 

a symmetry which causes derivations in the literature to differ in appearance. 
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The integrated radiance from all the volume elements along a given 
line of sight is from equations (6) and (16) 


* fo /a maX V x ” (a ' r > i dAda < 23 > 

mm 


where a and a are the minimum and maximum values of the particle 
min max 

size. 


Equation (23) has been derived in various forms by numerous 
authors; e. g. ; Searle [ 1] , Schoenberg [ 2] , Van de Hulst [ 13] , Walter [ 15] , 
Giese [ 26] , Gillett [ 27] , and Aller et al. [ 22] . 

C. Particle Size Distribution 

Due to the mathematical difficulties involved in solving Maxwell's 
equation for a scatterer of arbitrary shape, the approaches to the scattering 
functions are limited to using the following: 

1. Mie theory formulation for spherical particles. 

2. Formulations for other basic shapes, e.g. , cylinders and spheroids, 
from which the scattered radiation can be described in a closed 
form. 

3. Perturbation methods and approximation. 

4. Laboratory results. 

Considerable analysis is based on the assumption that the interplanetary 
dust particles are spherical, even though in actuality, spherical particles in 
interplanetary space are thought to be rare. On the other hand, a large number 
of spheroids have been found in recent Apollo finds. These spheroids were 
formed by a process which may have produced spheroids in interplanetary space. 
However, one cannot, a priori, extend the properties of particles in the finds 
to particles in interplanetary space, even though the two may have had similar 
origins and histories. 

In support of the assumption of spherical interplanetary dust particles, 
one must also note that distributions of large volume, irregularly-shaped par- 
ticles can be approximated by equivalent distributions of spherical particles [28]. 
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The distribution function r) (a, r) can be separated into a particle-size 
distribution n(a) and a spatial distribution s(r) .' The size distribution function 
n is arbitrary, but in zodiacal light model calculations the assumption of an 
inverse power law is dominant because. of meteor studies, physical considera- 
tions, photometric analysis, and simplicity. 

The inverse power law particle-size distribution (Category IV) 



will be primarily used here. Section II. D lists other distributions which could 
be used. However, the particle sizes and materials have to be chosen so that 
the observed -radiance and polarization curves are reproduced.. It may be pos- 
sible by using the dynamical forces (see appendix) and chemical and physical 
processes working on the interplanetary matter to derive a particular distribu- 
tion, although this has not been done. The forces and processes are important 
in the selection of a distribution form. 

D. Heliocentric Distributions of Scattering/Particles 

The spatial distribution function s(r) is primarily determined by the 
solar gravitation and radiation fields and, therefore, is assumed to be a= 
heliocentric distribution and symmetric about the ecliptic. The various 
heliocentric spatial distribution models are divided into planar and ellipsoidal 
models. The planar models, due to Gillett, are the constant thickness plane 
model and the linearly increasing cone model. These models represent first 
approximations for calculations of the photometer functions out of the ecliptic. 
A somewhat more realistic model, the exponential model due to Aller et al. 
[22], has an exponential decay factor for the spatial density of the particles. , 
The decay is in a direction perpendicular to the ecliptic. 

The other models have ellipsoidal symmetry and are, at least partially, 
based on the orbital motions of meteors. The flatness of the ellipsoidal sur- 
face is characterized by one of the parameters. 
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, Before studying these models, consider the photometric functions as 
observed^ in the ecliptic plane fop heliocentric models. ' 

Let the function r) (a, r) have the form: 

^■ r '= n o(r)' P (S) m («) 

where n (a/a ) _p is the particle size distribution and a , p, m, and n are 
o o o o 

parameters for the model. 

The radiant flux surface density has an inverse square dependency on 
-the solar distance due to the geometric attenuation: If we define the spectral 

radiant flux surface density at R =' l AU to be E°, then the general spectral 
radiant flux surface density is 



Equation ( 23 ) then becomes . - s . 



From Figure 11 (Fig. 10 modified) , it is seen that ' " ' - • - '■ ■ 

A = R cos r cos 9 , . .. ' ' " ' (27) 


where e is the elongation and © is the scattering angle. From the law of 
sines, we have 

r _ sin € 

R - sin © 


( 28 ) 
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hence, 


A = R cos € - R - S r— " cos 0 . (29) 

sin © v ' 

Therefore, 

dA = R sin edG/sin 2 © . (30) 


Substitution of equations (24), (28), and (30) into expression (26) gives 




Note that we now have an integration over © and that all scattering 
angles from e to ir contribute, as illustrated in Figure 12. 
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i v>\ ran be written ns 
Equation ( 32) can u 


BE® 6 vt N q 


sin © 


(— ) 

Van/ 


Van/ 


a max . / a \" p /bcrX 

- IT / n ° V%) W 

o a 


N = J 


a min 


n o (k) 


r .rt.de, per cubic meter at 1 AO and (*,) “ 

N is the total number of P _ over the distribution of particles. A 

Je averse diiierential cross secuo ^ **. 

ehange of the integration vartabie from 


— u a 

(*l> — ^ 

Van/ 


V an/ 


(33) for the simple behavior of Iy assume tha 
To investigate equation ( ) 

— f Hence, consiaer the integral, 

<£. is inaepenaent of 0 • Henc 
An 
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I 

m 


(O 


. m+1 . . 
sin (e ) 


7T 

r m . . 

J sin (©) d© 


For m = 0, 


I 

m 


■ _ 7 r - € 

0 - s in e “ 5 


for m = 1, 


I 

m = 


1 


1 + COS € 

sin 2 € 


and for m = 2, 



b (n-€ + sin e cos e) 
sin 3 € 


For m = .0, there is no radial dependence, (R/r)°= 1, except for geometric 
attenuation of solar radiation. These results are illustrated in Figure 13. 
Notice that for e = 45 deg, I is approximately independent of m . Also, 

I has the form of the observed radiance data which indicates in the ecliptic 

m _ dcr 

that the dependence of L on — is limited. 


As e approaches n , 


ii" , y 

. m , 1 

, , sm © d© 

. m+i -'e e .. _ . _ 

sin e y ° 


1 — f 

m+l Je 


— -H- - -hn • 

-~ir m+l •’ m + l 
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ELONGATION (DEGREES) 


Figure 13. I dependence on e for 
m 

m = 0, 1, and 2. 
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hence, I 

m 


( € = 7T ) 


— . This result can give an estimate for the value 

m + 1 6 


of the gegenschein, assuming that 
behaving smoothly in this region. 



is known and is 


E. Gillett Models 

The two models now considered which are not limited to the ecliptic 
plane are due to Gillett [27 ]: 

1. Model I - Constant thickness model (Fig. 14). The particles are 

distributed in space-like s(r).= (R/r) m up to a perpendicular distance W/2 
from the ecliptic plane. For perpendicular distances greater than W/2, 


2. Model II - Linear thickness model ( Fig. 15) . The particles have 

a spatial distribution of (R/r) m up to a perpendicular distance from the 
ecliptic equal to r sin and s(r) = 0 for distances greater than r sin 


Since the physical density would not have the discontinuity from 

s = (R/r) m to s = 0 at the boundaries of Models I and n, it would be desirable 
to use a combination of layers of width W lf W 2 , W 3 , etc. , for Model I and 
various cone angles cj etc. , for Model II. In addition, the components 

would not all have the same weighing factors, so that (e , i) would have the 
form: '■ 


- N 

L x ( e » o = Z 

x £ i 


k if 

A A 


(e» i) » 


(37) 


with the sum over the various layers. The angle i is the inclination of the 
scattering plane from the ecliptic. 

Figure 16 shows the parameters of Model I. At point P, where the line 
of sight is at a perpendicular distance W/2 above the plane of the ecliptic, the 
scattering angle 0 = 0 max « The length d, the distance from P to the earth- 

sun line, is given by d = W/[ 2 sin ( i)l . 0 , ib, e , and d are all in the 

’ max 

scattering plane defined by the lines of sight to the point P from the earth and 
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Figure 15. Linear thickness model. 




Figure 16. Geometry for constant thickness model 





the sun. The scattering plane is inclined an angle i to the ecliptic. In the 
scattering plane , ; 


0 = ib'+€ 

max 


( 38 ) 


and 


sin ip = d / \l d 2 + (R - d cot e) 2 ; 


(39) 


hence, for Model I, 


0 ( € » i> W) = €’ + sin' 

max 


1 (l/l 1 . + (— sin i- cot c) 2 (40) 


For Model II ( Fig. 17) P is the point which/is a distance r sin oj^ 
from the ecliptic plane. The s^1^r'iSg'aii^le"fdr;t))l8 point is e max ( e » U « ) 
The distance d of P from the earth-sun line is ^ 


d = r sin a; /sin i 
o 


/ £> 
* o' 

J V. 


v ' r f i 


(41) 


In the plane defined by the line of sight and the earth-sun line 


0 = 5 + to 

max 


(42) 


and 


d : /. . . . 

sin a; = ~ := sinu) /sin i ; 
r - o 


(43) 


46 




therefore, 


0 (c, i, w)=e+ sin 1 ( sin co / sin i) . ( 44) 

max o o 


From this equation it is seen that for i < co q , the line of sight never passes out 

of the dust envelope and no 0 < tt exists. For these values of i, 0 is 

max max 

taken as 0 = tt. The line of sight can be such that 0 is doubly defined 

max ° max ^ 

as shown in Figure 18 or as seen from the multivalueness of the sin in 

equation (44). In this case, the integration limits are (e , © ) and 

max 

(0* , vr) . The formula can be transformed to 

v max ' 


0 (e,, i, co ) = e + tan /sin to / \l sin 2 i - sin 2 co ) 

max o \ o o / 


for sin 2 a; /sin 2 i < < 1 . The double cross relation is 
o 


0' (e , i, co ) = e - tan 
max ' o 


-l 


sin u> 


nj { _ cin ^ 


sm‘ i - sin‘ u> 


+ TT 


(45) 


for e < 0’ <7 r. For i < co , then 0’ = tt. 

max o max 

The distribution of dust away from the plane of the ecliptic and up to a 

scattering angle, 0 » has the same form as in-the-ecliptic plane for 

max 

Models I and II. 0 is determined by the elongation e and inclination i. 
max 

The expression for the spectral radiance out of the ecliptic can be 
written as [ 27 ] 
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Figure 18. The double cross of the linear thickness 
perpendicular to the earth-sun line, the contribui 

e to 0 and 9 1 to it. 


RE A 6A n 

V 6 - ‘>° — 557 f 

sin £ e 


9 ( e , i) 

max v 


slnm0 (da) de (46) 


Define 


RE 8A N ‘ 
A o_ 

. m+1 
sin e 


'f sinm ® (da) de 


-/ 


. m 
sin e 


0 


max 


(?) 

d0 

V dfi / 



(47) 


\ (3R, O) 


RE 6A N 
A o 

. m+1 
sm 311 


7T 


/ sinme (55) de 

911 v 


(48) 


Now, equation (47) can be written in the form of in-the-ecliptic integrals: 


l a c*. i) = l, (« 


/ sin 0 \m+l 

, o) - ( — - g - ^ ) L,(9 . , o) 

' \ sm e / ' 


A max 


Consider the spectral radiance equation for the ecliptic plane (i = 0): 


7T 


Vo = \ °) = 


/m 

sin © d0 ( 50) 


m+1 
sin e f 



e 

where C = R E 5A N . Then, because 
\ o 

O < £ £ 7T, 
e < 0 ^ 7r, 
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and 


o < sin.0 '< 1, 


0 

\ dfi/ 

for any value of m> 0, a smooth which reproduces (e ) can be 


found as long as L ( e ) , sin 
A 


m+l 


€ generally decreases with increasing e . The 

integrand is positive and, as £ — ,it, the area under the curve decreases as 
e increases toward its upper limit. Thus, there is a wide range possible for 


manddu/dfl which will reproduce L (e). If, however, it is required that 
A 

dcr/dfi reproduce L (e , i) for points off the ecliptic, then another set: of 

A 

parameters, 0 (e , i ) , which describes the variation of the density of 

IHclX 

particles away from the ecliptic is introduced and the possibility exists of 


further limiting the range of acceptable values for m and d cr/dQ : 

0 U, i) 

( £ . ■) = — — — / 


L v 


max 


; m+l . 
' "• sin e 


. m 
sin 0 


Vdfi/ 


d0 


(51) 


Fo r the s a me valu e of m, an average polar izatiori function g (0) where g (0) 

= i 1 (0)-i 2 (0) can be found which will reproduce the integrated polarization 
in the ecliptic 


7T 




m 


„ _ .m + l 
2 L sin e 
A 


. J . sin 1X1 0 g (0 ) d© 
• e 


(52) 


This can be done because g(0 ) can assume negative as well as positive values. 

A 


This g(© ) together with © max ( € , i) must predict values of (e , i) : 
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— . m+1 

2 L ( e , i) sin e 
A 


r 


e 


max . m 


sin 0 g (0) d9 , 
(53) 


which are in accord with the observed values. 

The model calculations were outlined by Gillett [ 27 ] in the following 

fashion: 

1. For a given set of m’s and using 


\ ( e > o) = 


. m+1 
sin e 


r . m / dcr \ 

J sm 0 (STj 


de 


(54) 


(da/dfi) is found which reproduces the observed values of L (e , o) 

■ A • 


2. Using the (do/dfi) determined above, the values of 9 (e , i) 

mRx 

are determined by fixing 


U (« , i) 
A o 

L (e , 0 ) 
A o 



(55) 


where the experimental values appear on the left side of the equation and 
is a specific value. 


3. The values of m, (dcr/dfl) and 0 (e , i) are tested by determining 

lYlcLX O 

whether they can reproduce the values of the P ( e , o) and P ( e , i) . 

A A 

Gillett has used the difference 
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( 56 ) 


L (e, i) - L (90, 90) 

A A 

to compare observations of the zodiacal light. This allows observational values 
used to be less dependent on the total observed radiance and instrument calibra- 
tions obtained from various observers. 

F. Exponential Model [22] 

To avoid the discontinuities in Gillett’s models, the spatial distribution 
can be written as a smooth function. The exponential model assumes the form: 

s (2, Z) = N q exp (-KZ/R) 


where N is the particle density at the earth, v and l/K are density distribution 
o 

parameters, and 2 is the length of the position vector projected onto the 
ecliptic plane. 

( 

Since the geometry is the same as the ellipsoidal model which is to be 
discussed, Figure 19 can be used to define. 


r = R sin e/sin 0 , 

A = R sin (0 - e ) / sin 0 , 

, Rsine 

dA = : ~ 2 • ~ d0 

sin' 0 

Z = A sin 6 , 


and cos e = cos 6 cos f. . 6 and l are the ecliptic coordinates of the line of 
sight. If y is defined by 


Z = r sin y 
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then, from the equations above, 


R sin (9 - e)' sin 6 1 

sin 0 R sin e/sin 0 


and Z/r takes the form 


sin y = sin (9 7 e ) sin 5 /sin e . 


Again, from Figure 19, 

S/r = r cos Y/R = sin e cos Y/sin 0 


Therefore, 

1 

1 




With these expressions, the equation for s (S, Z) becomes 


s (2, Z) = s (e, 0, 6) 


sin 2 0 

2 

-Ksin (0 - e)sin 6 

sin 2 e - sin 2 (0 - e) sin 2 6 

exp 

sin 0 


(57) 


G. Belt Model [8, 30] 

In this model, interplanetary space ( in the ecliptic) is separated into a 
system of belts such that the boundary of each belt coincides with the orbit of 
a planet. The kth belt is populated with particles having a given size distribu- 
tion and any refractive index In each belt one must take into account the dis- 
persion forces which act or have acted on the particles. Behr and Sipdentopf [8] 
used the belts to determine electron densities versus distance from the sun. 

From equation (28) and Figure 21, we have 


0 . 
I 


. -l 
sin 


R sin e 

Y 


(58) 


so that the belt model components of the spectral radiance from equation ( 32) 
become: 


_ K 

h < £ - *> “ 


RE® 6X 

A 


. m+l 
sin e 


N 


l 

j=i 


q 

max 

/ 

a 

mm 



o 


x 



0. , 

r l+l . m „ 
J J sin 0 

0 . 

J 



d0 ( 59) 


where 0 = e. 

3=i 


Figure 20 indicates that equation (58) must be used with care 
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since a single can define two 9^’s. This model has the property that the 

spatial density can be increased at specified distances, such as at the asteroid 
belt [ 31 ] . 

H. Ellipsoidal Symmetry Model [32, 33] 

In the ellipsoidal symmetry model, the particle number spatial density 
distribution is assumed to have equal values at ellipsoid surfaces. These 
ellipsoids are in a Cartesian coordinate system centered on the sun with x- and 
y-axes in the ecliptic plane and the +z-axis pointing into the Northern Hemi- 
sphere (Fig. 22). 

For this model, the spatial number density s (x, y,z) of each particle 
component is assumed to decrease with the distance from the sun according to 


s (x,y,z) 



J 


(60) 


where in the ecliptic x 2 + y 2 =. p 2 , R is set equal to 1 AU and n^ is the number 

density of the particles at 1 AU. Outside the ecliptic, the particles are assumed 
to have an equal number of densities along the surfaces of the ellipsoids defined 
by 





(61) 


Out of the ecliptic plane, y is a parameter characterizing the flatness of the 

ellipsoidal surface corresponding to the number density, n^ , Wt>f [33]. In 

the ecliptic plane only, p represents the distance of the corresponding equiden- 
sity surface from the sun. 

A spherical distribution of the particle number density would correspond 
to y = 0, while values of y » 1 correspond to a strong concentration of inter- 
planetary dust toward the ecliptic plane (a disc) . 
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Using Figure 19 and equation (61), we have 
x 2 + y 2 + z 2 (y 2 + 1) = p 2 



where 


r 2 _ x 2 + y 2 + z 2 
Therefore, 



From the law of sines, we have 
r = R sin e / sin © 


and 


A = R sin (0 - e)/sin 0 , 

and the element of distance dA is given by 
dA = R sin e d©/ sin 2 © 


( 62 ) 


(63) 


(64) 


(65) 
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From Figure 19, 


z = A sin 6 (66) 

and 


cos e = cos 5 cos St . (67) 

Since we have a continuous distribution, no G exists; hence, the formulation 

max ’ 

of the spectral radiance equation becomes 


RE, SAN 7r 

L = — i ^ f 
\ (3+1 J 

sin e e 


R 

sin 0 


(— ) 

\dfi/ 


de 


'i + (Xi 


sin 6 sin (Q - e 
sin e 




P/2 


( 68 ) 


which is evaluated by numerical integration. 

As reported by Schoenberg [ 2 ] , Van Rhijn [ 3 ] used spheroids as con- 
stant density surfaces. He aLssumed the density distribution 


s (x, y) = n^ (R 2 - a (x 2 + y 2 ) - bz 2 ) 


with the condition 


R 2 -a(x 2 +y 2 ) - bz 2 = 0 


at the upper boundary of the radiance integral, x,' y, and z are the orthogonal 
heliocentric coordinates (i.e. , centered at the sun with the x-y plane in the 
ecliptic) . The integral cut-off condition can be written as 
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sin 2 0 - a sin 2 e (b - a) sin 2 5 sin 2 (0 - e) = 0 , 

max v v max ’ 


where 5 is the geocentric latitude of the line of sight. The last condition used 


x 2 + v 2 = R 2 sin 2 e - sin 2 6 sin 2 (0 - e) 

y .. L sin 2 0 J: 


and z as is in equation ( 66) . 

To determine the polarization at the celestialpoles, Fesenkov [ 34 ] 
assumed surfaces of constant density that were ellipsoids of revolution 


x 2 + y 2 ' " z 2 


Assuming that the size and spatial distribution are separable, 


V (a , r ) = n Q n (a) S.(r) 


and using the relations 


r 2 = R 2 +z 2 = R 2 ( 1 + g 2 q 2 - g 2 ) 


A = g n/ a 2 - R 2 = g 


R vq 2 - 1 , 


where g is a constant, q= a/R, and A has been assumed to be equal to z (i.e., 
along the line of sight), equation (26) becomes 


] E e A 2 <5AR 2 N ® s(r) i ; dA 
L A = 4? J r 2 J=1 ’ 
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As before, 
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A r 

N C - ST J 


max 
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n n (a) da 
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min 
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max 
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max 
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Q! . 

mm 


n (a) E (©, a) da 


f n (a) da 


a . 
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7 ) 


In terms of q, becomes 


A 2 E A 6 A R 2 N oo s (r) i (©) g Rq dq 

— / L 


h} = ^ 


47T 


\ R 2 (l-g + q 2 g 2 ) (q 2 -!) 1 / 2 


If s(r) varies as l/a, then 


L J = gC T N f 7 - 

1 ” n J M 4- 


L ( 0 ) dq 


o ^ (1 + g 2 q 2 - g 2 ) (q 2 - 1 ) I / 2 


where C f = E A 2 6 A R/ 47 T 2 . Assuming q= sec <j>' , then 
A 


gC' N 

7J = 21 

A R 


L 


7r/2 


i (0 ) d <p' 
•* 


2 2 

cos <f > ’ ( 1 + g tan $>’) 


where 


0 = cos 


' ( 


-g tan 0 
\/I + g 2 tan 2 


0’ / 



For any point in the sky, Fesenkov gives the general expression for the 
variables a, r, and 0 as: 

a 2 = R 2 + A 2 cos 2 5 - 2RA cos5 cosi + A 2 /g 2 sin 2 5 , 

r 2 = A 2 + R 2 - 2R A cos e , 
sin (0 -e ) = A/r sine, 

and 

cos e = cos 6 cos I , 


where 6 and l are the ecliptic coordinates for the line of sight. 

I . Meteor Model [35] 

As a logical starting point, observations of asteroids, meteors, and 
comets have been used to develop zodiacal light models. Most of the early 
work used large interplanetary matter as a source of the scattering functions 
or the equivalent phase function (e. g. , see Ref. 1) . In the meteor model to 
be described, the distribution of the orbital elements is used. 

In the meteor model, the concentration of dust material in interplane- 
tary space is calculated as a function of the distance from the sun and the 
heliocentric latitude $ on the basis of a bivariate distribution f(j, a) of 
meteor orbits with respect to inclination j from the ecliptic and semimajor 
axis a. 

i 

The number N( j, a) of orbits with parameters j and a within the inter- 
vals (j, j + dj) and (a, a + da) is 


N (j, a) = f (j, aj dj da . 


If the orbits are circular, then within a segment df of the orbit (Fig. 24), con- 
fined to the interval (4>, <t> + d$) in heliocentric latitude, a particle will make 
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X 


Figure 23. Parameters of an elliptical orbit. 
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TABLE 6. ELLIPTICAL ORBITS 
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a contribution to the number of particles in df which is proportional to the ratio' 
of the time dt it is within this interval to the total period T. This contribution is 
equal to the ratio of the element df to the circumference 2wa; i e. , df/27ra = 
dt/T. The element d£ can be found in terms of j using (Fig. 24) : 

I 

sin j __ sin (tt/2) 
sin $ 7 sin (A. /a) 

and ' . 

® = j? /a ; , t 


hence, 

sin $ = sin j • sin (i/a) . 3 • (70) 

Differentiating both sides with respect to f , we have 


d sin $ d$ . d sin l /a 


d£_ 

a 


sin j n/ 1 - sin 2 4> df 

sin 2 j a 


, d$ 

cos <f> — = 
df 


d$ df J df 


cos 


, d$ ■ £ X 

— = sin j [ cos — ) 
df J V a / 


By using equation ( 70) , we have 
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d* 

a 


cos # d# 


( 71 ) 



Assuming that the number of particles is proportional to the number of 
orbits and that the longitudes of the ascending nodes are distributed uniformly, 
then the number of particles located in the segment df and having orbits with 
parameters j and a is proportional to the quantity 


LLk*I 

2na 


dj da df f (j, a) cos # dj da d# 
2 7r \j sin 2 j - sin 2 $ 


where all these particles refer to a volume of dv = 2vr 2 cos# dr d# . The 
volume element dv is determined by the surfaces between #, $ + d#, r, and 
r + dr. A unit volume at a point having the coordinates (r, #) will contain a 
number of particles proportional to the quantity 


cos # f (j, r) dj dr d# 
2n \[ sin 2 j - sin 2 # 


1 

27rr 2 cos # d# dr 


f(j. r) dj 

47 r 2 r 2 \f sdr^'T^ldn 2 '# 



(72) 


Integrating with respect to j, we obtain the dependence F (# , r) of the particle 
concentration on the heliocentric latitude # and the radius (r = a) o'f the cir- 
cular orbit: 


/ 

# 


7 r - # 


j (i. r) 


's/sin 2 j - sin 2 # 


dj > 


F (# , r) = 


4 7r 2 r 2 


(73) 



where C is the normalization constant. Recalling equation (44) , 


0 = e + sin 


-l sin <£ 
sin i 


one finds that 

, #(0, e, i) = sin [sin (0 - e) sini] . . ' . 

From equations ( 28) and ( 30) , the integrated radiance [ equation ( 26) ] becomes 


CE X 6\ , max 


it - sin - ' (sin (0 - c) sin '1 


0J-. .. 

t (1, R sin e/sln 8) dl da <39 


\ (e. o = 4ir 4” s i'n T 7 /* fa . rJ,- -x ... n(a) sin J 0 dn ^ ^ -~ r - 


a min sin" 1 [sin (0 - c) sin i] 


nT sin 2 j - sin 2 (0 - e) sin 2 i 


da dO (74) 


where 




n (a, r) = n (a, j, r) , 


and n(a, j, r) is assumed separable into a particle size and a spatial distri- 
bution; i. e. , 


a (a, j, r) = n (a) n (j, r) 


This method can be generalized [21j by writing the distribution f (j, r) 
as a function of the following parameters: 

. I . . ' . ' . 'i 1 

a = semimajor axis 
e = eccentricity 

'•'« j = inclination ; ! - • ... 

s = radius of dust particle ....... 
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B = shape parameter 
5 = physical density 

and time. Then 


£ (q, e, j, s. B. 6. t) = d 6 r) (t)/da de dj ds dB d6 
1 1 J 


where t] refers to the total number of particles in interplanetary space. For the 
spatial and size dependency, it could be assumed that f is separable into a pro- 
duct of functions of each independent variable: 


f (a, e, j, s) = f t (a) f 2 (e) f 3 (j) f 4 (s) 


Assuming f 4 = s P , then the number density as a function of r and 9 is 


, CO 1 ' • 7r/2 

n ( r ’ 0 ) = ~ I I I 

r/2 I r/a - 1 1 tt/2-0 


~2 

a f(a,e, j) da de dj 


H-)'H 


sin 2 j - cos 2 0 


U/2 


(75) 


J. Generalized Method For Meteor Models [36] 

The particle number density can be generalized and formulated in terms 
of the distribution of orbital elements semimajor axis, a, eccentricity, e, 
inclination, j, node, 0, and argument of perihelion, co. The number density of 
particles having orbital elements between (a, e, j, £2, co) and (a + da, e + de, 
j + dj, + dfi, co + dco) with period P and velocity V is equal to 


d*N = 


f (a, e, i, Q , co) da de dj df2 dco 
r 2 dr sin© d0 d <p 


dL 

VP 


(76) 


The time spent by a particle in the volume is equal to dL/V where dL is the 
pathlength of the stream in the volume element at r, V is the stream velocity, 
and P is the normalization factor. 
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The differential area dft dco is related to the differential area d9 d</> by 
the Jacobian: 


dft cko 
dG d <p 


J ( ft , CO , 0, 0) = 


/ 9ft\ 
\9©A 

(&) 

/9ft v 

( — ) 

u; e 


(77) 


To deduce the Jacobian, refer to Figures 23 and 25 where the geometry of an 
elliptical orbit is defined and to Table 6 which gives the properties of an ellipti- 
cal orbit. The elliptical orbit is defined by 


a ( 1 - e 2 ) 
1 + e cos £ 


(78) 


where angles referenced in. Figure 25 are 

£ = true anomaly 

ft = longitude of the ascending node 

co = argument of perihelion 

r, 0, <{> = spherical coordinates 

i = angle of inclination 

$ = co + £ = argument of latitude 

co = co + ft = longitude (amplitude) of perihelion 
ft + co + £ = free longitude. 

Using the following basic equations of spherical geometry (Fig. 26), 
cos a = cos b cos c + sin b sin c cos A 

sin A _ sin B _ sin C 
sin a sin b sin c 

sin a cos B = cos b sin c - sin b cos c cos A , 
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Figure 25. Generalized meteor model geometry. 






Figure 26. Definition of the spherical angles a, b, c and A, B, C. 

J 

we have the following relations: 

cos (a; + f ) = sine cos (0 -O) , (79) 

cos 0 = sin j sin (co + g) , • (80) 

sin (u> + f ) cos j = sin 0 sin (<p ~ U) ( 81) 

and equation (81) divided by equation ( 80) yields 

cot j = tan 0 sin (<p - Q) . , (82) 

From equations (78) and (79), 

co = cos 1 [sin 0 cos (<p - S2)] - cos 
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but equation (82) yields ' •- 


.• -l 


fi = -sin (cot j cot 9) + 0 ; 

4 * t j • 1 

- * 1 * 

thus, the expression for co becomes 


co = cos 1 |,s in 0 cos [sin -1 (cot j cot 0)jj , - cos -1 “ gj 


-l 


* . o , ox 1 / 2 .-i Taf 1 - e 2 l 1 1 , 

co = cos (sin 2 0 - cot 2 j cos 2 0) -cos *• ,.> ■■■ ' ; (83) 

re e 


hence 


&}) .'o' 

Sf) e 


m . i 

8 ^e 


rc ■ 


'■ o 


0 


9£2\ _ cot j/sin 2 © 1 

99 / , ~ ( 1 - cot^j cot 2 ©) 1/2 

0 


( sin 2 0 - cot 2 j cos 2 9) ^ 2 (2 sin 9 cos 9 + 2 cot 2 j cos 9 sin 9 ) 
~ [ 1 - ( sin 2 0 - cot 2 j cos 2 0)] 1/2 


-sin 
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This dependence is seen directly from Figure 25. Therefore, 


J (O, w, 0, 0 ) 


sin Q/\j sin^j - cos 2 0 


( 84 ) 


The relations 


~ cos 1 u = -1/ ( 1 - u 2 ) 1//2 
du ' ' 

arc sin x = ± arc cos ( 1 - x 2 ) 


with a, e, j, and r assumed constant were used in the derivation. 

Note that dL/dt = V , hence, upon dividing the numerator and denominator 
of equation (76) by dt, one obtains: 

_ f (a, e, j, SI, u) dade dj 
r 2 P(a) (sin 2 ] - cos 2 0) 1/2 

From the dynamics of two particle orbits 


P = 27ra 3/2 /(GM ) 1/2 
* o 

and 


V 

r 


G M 

o 


e 2 a 2 - (r - a) 2 
r 2 a 
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Therefore, 


N = 


-i- r 

2lrr r/2 


1 7r/2 

/ / 

Ir/a-ll ir/ 2-0 


f (a, e, j) da de di 


a 


T 2 /r \ 21, . 2. 2 J 

e - lj (sin j - cos 9)< 


il/2 


(85) 


The argument of perihelion and the node is assumed to be uniformly distributed 
over(0, 27r). Azimuthal symmetry exists in dN(r), since 0 does not appear 
in the equation. 

If only gravity and radiation drag are important, the inclination of the 
orbit does not change with' 'time' and one can assume that • * .. . . 

f (a, e, j) = £ (a, e) £ (j) . , 


The integration over the inclination can be performed and will yield a quantity 
M(0), depending only on 0 : ' 


7r/ 2 

M(Q) = / 

7T/2-0 


L.(j)--dj ... . 

, . 2 . 2 ' • . 1/2 
(sin j - cos 0) 


( 86 ) •, 


If f ( a > e) = g(a) h(e), where h(e) is continuous between 0 and 1, arid if g(a) 
is given by a power law [If £ (a, e) = ea 2 for small e, then dN is independent of 
r]. 


g (a) ^ 


(87) 


then, upon making the change of variables x = r/a, 


dN 



h (e) de 

"2 , . 2 1 1/2 
_ e - (x - 1) 
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thus yielding a power law in r. 

For the particle radius s, assume 


.. '• d N (r) -••• ; o n / 2 ■ 

; - , c fc h w rV 


iffliL 


2 . 2 1/2 
tt/ 2-0 (sin j - cos 0) 


( 88 ) 


‘ .y : ' . ’ 


where 

!'<}■ ■■■ ■ 


f (a, e, j) = const, s P a Y h (e) f (j) r , 
s 


(89) 


2 - y 



and C , 
£,h 


(/ 3 ) is a constant depending on j and h. 


. ‘ ? £ (j) represents the distribution of inclinations on a sphere, but it has 

already been assumed that there is no <p dependence. Hence, a uniform distri- 
bution of inclinations would not give a uniform dust population on the celestial 
sphere. If the intersection of the polo of the Orbit ! With the sphere is represented 
by a point, then the density of points is given by 


# points _ # points 

dA - sin j dj d£ 


But, since there is no <p dependence, £ (j) must be proportional to sin j to have 
a uniform distribution of dust. 

, . ■ Distributions of the forms ,£ (j) =* exp(-bj) sin j 

£ 0) - (cos j) sin j j v 

(k = integer) are suggested from meteor studies [ 36 J . 
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SECTION IV. INVERSE FUNCTION ANALYSIS 
A. Integral Inversion 


The analysis of the observations of radiance L and polarization P of 

' A. - ' A, 

zodiacal light would be simplified if it could be reduced to a consideration of the 
average scattering intensity functions, where the average is over the particle 
size distribution. Several methods for doing this are considered in this section. 

Assuming (a) the particle size distribution to be independent of the 
spatial distribution, (b) one refractive index for all the particles, and (c) the 
spatial distribution to vary as (R/r) m , then the components of radiance in 
the ecliptic plane are given by 

= ■■ ■ ■■ nT - H ’ " " sin m 6 4(0) d0 
snr u e J € 


l a - r sin m 9 i 2 (0) dO’ 


where iT(©) are the average scattering intensity functions and C is* a constant. 

J v ’’ ; l -■ 

' .! ' i a .. - 

- v The radiance and polarization values are given by 






from which the relations between the calculated values (L 1 , I? ) and the ob- 

A A 

served values [F ' (e) , F<f (e)] can be .found: 

1 
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(90) 


L * = L X< 1+! V E V<«> 

= Z x (i -P X ) = F 2 '(€) . ’ (91) 


F ' and F 2 ' are the values obtained from observations and are related to 


l U 2 (0) b y 


F 


i» 


C 


. m+i 
sin e 


7T 

J sin m 6 ij, 2 ( 0 ) d 9 
e 


(92) 


For later convenience, define F ls2 (e) 


_ , , . m+i 

, X F U 9 U) sin e 
i> 2 ^ 


such that 

7T 

■ ■ ■ 

sin 9 i ± , 2 (0) d0 
e 


(93) 


Given Fj, 2 (e) from observations, the problem is to find i it 2 ( 6 ) . Equations 
(92) and (93) are of the form 

x 

F (x) = / K U, x) f U) d£ (94) 

a 


which are known as Volterra equations of the first kind, where K(£,x) is 
the kernel, i; is the auxiliary variable, and x is the current variable. Numer- 
ical methods for the solution of this type integral equations are given by 
Hildebrand [37] . Several such methods are presented below. 

B. Numerical Approximation Methods [37] 

A general method for obtaining solutions of integral equations consists of 
reducing the integral to a finite set of algebraic equations. The solution of 
equation (94) 
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F (x) = jk K (x,4) y ($) d| 


(95) 


may be approximated by a linear combination of n suitably chosen functions 

<p , <p , . . ., 0 of the form 
12 n 

n 

Y < x >“ E A k*M x ) . (96 

k=l 

The n constants of combination are determined to satisfy equation ( 95) as 
nearly as possible over the interval (a,b) . If y(a) and y(b) are known, it 
is convenient to assume the approximation in the form 

n 

Y(x) - (p (x) + Yj A k W 
k=l 

where <p is chosen so that it satisfies the end conditions, 
o 


( a) = y (a) 
^(b) = y (b) 


and the remaining (p's are made to vanish at the corresponding ends of the 
interval. In some cases, (p may be chosen so that it approximates y(x) 
closely. The remaining terms, the summation, can then be considered to 
constitute a perturbation expansion. 


C. The Method of Collocation [37] 


For the particular problem, equation (93), assume that 


n 



1 > 2 ^ a k’^ 


(97) 


where the cv^'s are preselected particle size parameters and the C^'s are 
constants whose values are to be determined. Then 


n 


F l» = Jc sinm (^ 2 c k 2 (a 9) do 

k K 
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becomes 

n 

P lf 2 (c) = Z C k sin m (0) i t> 2 (a k ,fl) d6 (98) 

k 

= Z c k 4’ 2l6) • (99) 

k 


A set of n equations for the determination of the C 's can be obtained by 
requiring that equation (99) be an equality at n distinct points in the inter- 
val (e, 7 r) . If we denote these points by e., the resultant equations are then 1 
of the form: 1 


X C k g£ 2 ( €i ) =F 1> 2 ( ei ) =F 1>2 (€.) i = 1, 2, 3, . . . , n 
k=l 


or in matrix form 
C G = F . 

D. The Methods of Least Squares [37] , ... . 

The method of least squares requires that the integral of the differ- 
ence between two functions over an interval be as small as possible. In the 
present case we wish to determine the best fit of 


f 17 sin m (0) 
6 

n 

Z C k g k (fl) 


k=i 


( 100 ) 


to the observed values F (e) by varying the C k 's for a chosen set of known 
functions g^ (0) . Hence the method of least squares requires 
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r 

J o 


n 

E C k /J sinin ( 0 ) Sk ( 0 ) de - F (e) 


k=l 


de 


(iOi) 


be a minimum. For notational simplicity, define 
<P. (e)= J n sin m (0) g k (0) d0 


( 102 ) 


In order that equation (101) be a minimum, its derivatives. with respect to 
each parameter C k must vanish; thus, ' ■ 


r.* t 

’n 

E c k^k( e ^- f <- e > 

o 1 

Jk=l 


or 


k-i 


' de = 0 i = , 2, ... , n ■/ (103) 


C k 1**1 <*) \ < e) d€ = fo* i (e) F ( e > dc • 1 '(104) 


Since F(e) are observed values, this is not a functional form, and a numerical 
approximation method must be used to solve equation (104) . The integrals 
are approximated by a weighted sum of the relevant ordinates of N conve- 
niently chosen points resulting in 


n 


N . ■. N • ; 

E C R E D r ^(e r ) <t> k (e T ) = E D r <M*r> F (e r ) 

k=l r=i r=i 


i = 1,2, . . . , n 


(105) 


where D (r=i,2..N) are the appropriate weighting functions associated with 
points €j, e 2 » • • •» involved in the numerical integration method. By 
suitably choosing g k (0) , the integrals of equation (102) can be carried out in 
closed form. 
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Several integration procedures applicable to equation (104) exist; e.g. , 
a Gauss-Legendre quadrature, the trapezoidal quadrature, Simpson's rule, 
etc. For the examples given here, Simpson's rule was used. Simpson's rule 
is based on the area under a parabola and is, 

j^y (x) dx = + 4y, + 2y 2 + . . . + 2y R _ 2 + 4y n _ 1 + y n ^ 


where 


yo = y (a) 

y t = y (a + h) 


y n = y + nh; ~ y (b) 


and h = (b-a)/(n-i) for n>2 and n odd. But, since equation (105) is homo- 
geneous in the D’s, the coefficients can be taken as 2, 1 1, 2, In 

matrix notation, equation (105) becomes 

a C = P (106) 

where 

N 

a ii = Ti 0ri D r <t> H i, j = 1, 2, . . . , n, 

J r=l 

/ 

N 

p. = T. <t> . D F , 

LJ ri r r * 
r=l , 

4> ri = 0i (€ r ), (109) 

r ( • ’ * 

and 


(107) 


(108) 


D r = D. 

r i 


6 .. 


ir 


( 110 ) 
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Equations (107) and (108) can be written as 

< . * ' ' 

<* = £ T D <$>_ ’ (111) 

£4 T DF (112) 

where, since D is diagonal, 

T T 

± D = (D£) . 

T 

T signifies that a matrix has been transposed, [e.g. , (<£ ..) = cp ..] . Hence, 

equation (106) using equations (111) and (112) is the matrix equation that has 

to be solved to find the C! s. 

k 

, When N = n, the number of chosen points is equal to the number of 

C. 's to be determined. However, when N > n, one can choose a number of 
k 

points greater than n and require that the integral equation be satisfied as 
nearly as possible at those points, thereby obtaining a best fit for all points. 
The errors committed by not satisfying the equation at the respective points 
are weighed in proportion to the influence of the ordinate in the integration! 
of the squared error over (0, 7r) by the coeffecients D. . The comparison of 

F (e r ) 

and 

n 

/ V sin m (0) £ C k g k (0) d0 
€ r k=l 

at the N chosen points gives an indication of the accuracy of the solution. 

The result gives i t , 2 (0) with the assumption of the spatial variation 
and number density; hence, uniqueness is dependent on these two parameters. 

For further analysis of i (0) = g^ (0) , one can use the least-squares 

method a second time so that 
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/ 


i (e) 



max 


a . 
min 


n 

Z 

k=l 


a. a 
k 


-k 


(0, a) da 


d0 = min 


where 





is an assumed form for the particle size distribution. Possible candidates for 
the values of various parameters such as refractive index, a , a . , and 
n can be found by this process. 1 

The methods of collocation and least squares were used for the inten- 
sity functions 

j 

Model A: i ^ 2 (0) = a^’ 2 + a 1 ’ 2 e -6 

Model B: l lt 2 (0) = bj’ 2 + bj’ 2 e~ 6 + b*’ 2 sin 9 

— 12 — 3 

assuming a particle density of N =10 cm . The results are given in 
Table 7. m in equation (98) is set to zero. In the collocation method, 
equalities were required at e =40 and 120 deg for Model A and e = 40, 120, 
and 160 deg for Model B. Weinberg's [17, 18] data were used as the experi- 
mental data. The intensity from Model B is the better fit to the data in both 
collocation and least squares than that from Model A. Both the collocation 
and least squares fit tend to be inaccurate for angles larger than 140 deg. This 
would indicate that data for these angles are not weighted enough in either 
method. 

Neither method gives satisfactory polarization for the models used. 

This is the main problem in using average scattering functions; i'.e., either 
the polarization or the intensity can be matched easily, but it is not as easy to 
find one function that will reproduce both. It is apparent that these functions 
are not unique . 
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E. Liouville Solution [37] 

Relaxing the assumption that the particle-size distribution is indepen- 
dent of the spatial distribution and allowing the particle- size distribution to be 
dependent on the spatial distribution, then another inverse method is available. 
By using 

dF(x) _ rB(x) 9K(g,x )f(£)d£ + K(x, B)f(B)~“ - K(x, A)f(A) ~ , 
dx “'A(x) 9x 


and if i^ 2 ( 0 ) = i u 2 (e ,6 ) , then 


9e 


-sin m e i (€=A, e) + f 1T si 

J € 


. m. 9 i ( e. 0) 
sin 6 — ~~ — d0 
9e 


This can be used only if i 1>2 is dependent on e, which will be the case if the 
particle-size distribution is spatially dependent; i.e., 


a 

„ max 

i-u 2 ( e > 8 ) = J v (e, cs) ij, 2 (a ,0) da 

a 

mxn 


Solving for i{ e, e) , 


-9 F ( c ) 

i(e,e) = 9c 
sin m € 

to which one can apply the method of successive substitutions; i.e., assume a 
value for i on the right side and solve for the left side, then iterate upon the 
solution . 


X 


sin^e 


7T sin 6 9i (e, 0) 

9e d0 
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F. Southworth Analysis [38] 


Southworth [38] analyzed, by a least squares method, the data of 
Smith, Roach, and Owens for the phase function and pa rticl e size distribution. 
It is assumed that a unique average scattering function i(fl) exists and the dis- 
tribution is heliocentric of the form 


77 (r,/3)= n(r) s (/3) , 


and symmetric with respect to an axis through the sun and ecliptic poles. /3 
is the heliocentric latitude and n(r) is defined by 

' log 1() n (r) = log 10 r + C 2 (log 10 r) . • (113) 


The phase function and the function for the latitude distribution, S, are taken 
to be the linear combination of seven triangular functions. 


and 


7 

s (sin/3) = Xs jqj (sin/3) 
j 

7 

Tie) = IX ‘k {e) 

k 


(114) 


(115) 


( 0 ) is defined by tk ( 0 < 0 k-l^ = °» \ ( 0 > 0 k+l^ = °» and fc k ( 0 = 0 k) = 1 . 
t (9) between ^ and e ^ and between 4^ and is represented by straight 



L 


7 

/ nsidf = H s 
line of j J 

sight 


Z 


k 


/ 

line of 
sight 


q. t, n d£ 
J k 


(116) 


Southworth stopped the integration at r =5 AU since contributions beyond this 

are negligible. A least squares successive approximation method was used to 

solve the equations for q. and h, . For Southworth's case with the least root 

J k 

mean square residual intensity, 
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log 10 n(r) = -1. 


Slog 10 r 


- 2. o( log 1Q r) 


The phase function and latitude distribution at selected values of sin j3 and 9 
for the triangular functions are given in Table 8. For 9 > 54 deg, Lambert' s: 
law roughly fits the phase functions. For 9 < 54 deg, a diffraction peak appears. 


TABLE 8. SELECTED VALUES OF i (0) AND s (sin 0) 


9 deg 

30.0 

54.0 

79.2 

104,4 

129. 6 

154. 8 

180. 0 

TJe) 

0.52 

0.00 

0.26 , 

0.40 

0.64 

0. 76 

1.00 

i » 

sin )3 

0.0 

0.1 

0.3 

0.5 

0.7 

0.9 

' 1.0 

s ( sin j3) 

1.00 

0.84 

i 0. 51 

0.38 1 

0.25 

0^ 20 

0. 23 


G. Cumulative Function [39] 


i 

Powell and Circle [ 39] formulated a mathematical procedure to remove 
the integrals over the scattering angle and particle size which occur in analyz- 
ing observations of the zodiacal light. The resulting algebraic equations are 
expressed in terms of scattering functions with known and recognizable prop- 
erties. Consequently, the uniqueness of deductions concerning the nature of 
the interplanetary dust is determined by studying the mathematical behavior 
of these functions. The development in this section is based on that given by 
Powell and Circle. 


Consider a uniform size distribution of particles N (X,R^) per unit 

volume per increment 6 a in the range 0<Q!<X at a radial distance R q from the 

sun. X is called the termination parameter since it terminates the unit step 
function at 


2ns. 

max 



(117) 


where a is the radius of the largest particle in the uniform distribution, 
max 

The scattering from the terminated unit step size distribution of particles with 
refractive index m* is given by the cumulative scattering ratio: 

X 

F.(0,m*X) = / o i. (9 , m* a) dor , 3 = 1,2 . ( 118) 
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The explicit functional dependence of m* will not be shown in the following 
equations. 


• t An arbitrary, continuous size distribution of particles, n (a) , can be 
expanded as a sum of terminated unit step functions with appropriate coeffi- 
cents, C (X ):' ■ 


a 


a (a ) £ C(X n ) N(X n ) , 


n 


(119) 


’.where the radial dependence of N(X, R ) has been assumed and not shown 

. o . . 

| explicitly. Some of the coefficients C( X^) may be negative, indicating a lack 

' of particles in the range defined by X^. However, the sum of the negative 

• coefficients must be less than the sum of the positive coefficients since a 
.negative number of particles in any region has no physical meaning. Using 
the step functions, the scattering from the actual size distribution of particles 
can be written as ... 


00 


L n (a) i (9, a) da . £ C(X ) F ( 6, X ) . , . . . , ( 120) 


•J 


n 


n j n' 


The spectral radiant equation for ,a polarization component j: 


_"j R.E®,6A A 2 ; 7T 

L A = m+ i' ■ ' / Sin m 0d0 f n(a) i (0, a) da 

87t z sin e e o i 


( 121 ) 


becomes 


_ j RE'- 6 A A. ir 

- L-, = ; — — — — T C(X ) f sin m 0F. (9,X ) d9 . ( 122) 

■ A * o . in+1 1 ■ u v n' * j v * n 7 v A 

8ir sin € n e . • 

The cumulative polarization function, p(0 , X ) is defined as 

m - . 

p,(0.x) - • .U*) 


J,;. 

where 


F(0,X) = F^e.X) + F 2 ( 0,X) 
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The measured polarization, P , due to all terminated step functions required 

A 

to characterize an arbitrary continuous distribution, is given by 


£ c(x n ,R o ) f sin m 0 A(0,X n ) d0 
; = -5 L 

X X C ( X n* R o ) f shim e F (e ' X n ) de 
n e 


(124) 


where 


A(0,X n ) = F 4 (e,X n ) - F 2 (0,X n ) - p(e,X n ) F(0,X n ) . (125) 
Hence, one has the following relation: 

fr- = Z c ( x n ) e A ( 0 * x n ) d0 

n e 

(126) 

Now define two functions of the measured radiance and polarization: 

M(m,e,A) = L (e) sin m+1 e ( 127) 

A 

and 

W( m,e , A. ) = P. (e) M(m,e,\) ' (128) 

A 


L * P M 


: m+1 
sm e 

RE 6 6 X 
A 


Then, 

M(m,e,A) = K(A)£c(X n ) /Jsin” 1 0 F ( 0,X n ) d0 

n 

and 

,’W(m,e,A) =K(A)^C(X n ) j\in m 0 A (Q,Xj d0 

n 


(129) 


(130) 
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where 


K(A) 


R 6 A A 2 
A. 

Stt 2 


Let 


7T 


J sin m 9 F(d) &6 = q ( tt) - q(e) 


and 


J ^ sin 6 A(0) d0 = Q(?r) - Q(e) 


Now, differentiating M and W with respect to e , and rioting that 


-dqUI = 9 Q(?) = 0 

de de 


one obtains 

9M 

de 

and 


= -K(\)£c(x n ) |a- 


n 


9W 

de 


= -k(a)£c(x) 


n 


9Q 
n' de 


However, 


9q . m ... 

= sin e F(e) 
9e ' 


and 


9 Q . m . . 

— = sin e A (e) 


( 131 ) 


( 132 ) 


( 133 ) 


( 134 ) 


( 135 ) 


( 136 ) 


( 137 ) 


( 138 ) 
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so that 


9M 

9e 


= - K( A) sin m e 2 c ( x n ) F ( e » X n ) 

n 


9W 


9e 


i = -K(a) sin m e £ C(X n ) A 
■ n 


(139) 

( 140 ) 


From equations ( 125) -‘and ( 126) , 


9M 9L A^ £ ^ . m+1 , v— /■ N . m 

— — = — sm e + ( 1 + m) L (e) sm e cos e (141) 

o€ 06 A 


and 


8P *0) . - 

37 = ~3T~ V £) SB £ + \ 


;9W 


9L. 


9e 


A . m+1 
sm e 


• •' • + (Ism) L sin e cos e 

....... . A 

•• • -v • .. <• V .. . + /• :v ;t • ■*<•' (.142) 

Now define two new functions of the measured radiance and polarization: 


‘ S(e,A,m) = 


: T(e,A,m) ■= 


-K 1 9M 


. m 9e 

sm e 


- K 


_1 


9W 

. m 9e 
sm e 


(143) 


(144) 


Combining equations ( 137) through ( 142) , one obtains 


S(£,C,X n ) = -2c(X n )F(e,X n ) 
n 


(145) 
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and 


T(e,C,X n ) - - Z c ( x n ) A ( £ « X n ) 
n 


( 146 ) 


S and T are obtained from the measured radiance L and polarization P for 

A A 

various assumed values of m as follows: 


S(e,\,m) = 


-1 


K(X) 


( 1 + m) cos e L (e) + 

A 


de 


sin e 


(147) 


T(e,X,m) = 


K(x) 


9P ' 

L x (e) ( sine ae 


_ 3L 

+ p x sr sta e 


+ ( 1 + m) P cos €| 

A 


(148) 


Equations ( 145) through (148) represent the fundamental mathematical rela- 
tionships between optically measurable quantities and the physical properties 
,f the interplanetary dust. S(e) and T(e) can be obtained directly from the 
observed curves and plotted for various values of m. The resulting curves 
are smooth and continuous. However, the measured curves, T(e) and S(e) , 


obtained by combining L , 

A 


d h - . 8P x 

— ■ p x • and 


, according to equations 


(147) and ( 148) must be matched by a sum of the theoretical curves, F(e,X n ) 

and A (e.Xj , as in equations ( 145) and (146) , with appropriate coefficients, 

C( X ) , to determine the unknown X , C( X ) and the refractive index which has 
n n n' 

been implicit in the equations. Powell and Circle accomplished this as follows: 
1. Characterize the measured curves S and T. 


2. Investigate the behavior of F(e) and A (e) by varying the refractive 
index and X . 


n 


3. Choose a refractive index and X for which the curves of F (e) and 


n 


A (e) are similar to S and T. 
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4. Substitute F(e) and A (e) into equations (145) and (146) and solve. 

5. Find a set of curves, A (e,X ) , which when summed match the 

» ' n 

product C(X ) A (e,X ) , corresponding to nonnegative coefficients, 
n n 

This scheme is a systematic way to match scattering diagrams. Steps 
2. and 3. involve a great amount of computation and study since, for a large 
spread in o', the curves are complex. Also, the fact still exists that this pro- 
cedure does not give unique values for m* or X . 
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SECTION V. OBSERVATIONAL RESULTS 


Although the radiance and polarization of the zodiacal light have been 
measured many times, there is considerable disparity among the values found 
by observers. Most observers have restricted their studies to the deduction 
of radiance of the zodiacal light in the ecliptic as a function of elongation from 
the sun. Observations at small elongation are necessarily close to the horizon 
and although this region is intrinsically the brightest, it is more subject to 
uncertain corrections for atmospheric extinction and scattering. Larger 
elongations can be observed at greater distances above the local horizon, but 
uncertain corrections for extraneous light then plague the reduction procedure. 
Most investigators regard it sufficient to assume that the zodiacal light is 
symmetric about the sun and the ecliptic plane and have reduced their measure- 
ments accordingly. According to Weinberg [ 40] , the main factors contributing 
to the differences in the reported values are as follows: 

1. The lack of dark-sky, low-latitude, high-altitude sites. 

2. The lack of a proven method for separating zodiacal light from air- 
glow and starlight, especially at high geographic latitudes and when using 
broad-band spectral systems. 

3. Difficulties associated with absolute calibration. 

4. Limited observational coverage, in time, over the sky. 

5. The lack of a satisfactory formulation for the effects of tropospheric 
scattering, especially near the horizon. 

6. Real changes in the zodiacal light, if any. 

The most photometric observations measure (l) the total radiance 

(L : ) , which for ground-based studies is a sum of various sources, (2) the 
obs ' 

total degree of polarization (P^) * which is dominated by the zodiacal light, 
polarization, and (3) the orientation angle of the plane of polarization. 

The measurement of the total radiance involves a number of terms 
from which the zodiacal light must be separated [ 17, 18] : 
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obs 


"^zodiacal + "^galactic + "^integrated 


L light 

+ L 


-n M t (z) 


scattered 


light starlight 

^ ^airglow + 


airglow 


components L emission lines continuum J 


>— t 2 M 2 (z) 


(149) 


! where t u t 2 are diffuse -source extinction coefficients for the astronomical I 
and the atmospheric components, respectively; 


M t (z) is the total air mass (equal to unity in the zenith at sea level) ; 


M 2 (z) is the air mass to the level of the atmospheric components; 


L . , is the component arising from the line and continuum emission 

airglow 

earth’s atmosphere; 


» L ^ is the component arising from scattering in 

scattered components 

the earth’ s atmosphere; 


L. , J is the diffuse component of radiance that comes 

integrated starlight 

from faint stars which are below the instrument threshold for resolution; - 


^galactic light is t ^ le scatterecl light from interstellar grains (which 
would be expected to have a polarization component) ; 

L .... is the light from the interplanetary particles, 

zodiacal light 

-Various methods are invoked to determine the zodiacal light from the non- 
linear multivariable function L f 41 1 . These methods usually involve the 

symmetry characteristic of each source to its natural coordinate system (e.g. , 
ecliptic, altazimuth, and galactic coordinates) . . 


Figure 27 is a collection of representative results for the zodiacal 
light radiance in units of S 10 (vis) between elongations of 25 and 110 degrees 
in the ecliptic. The data included are only in the visible spectral region. 


The relative variation of the. radiance of the zodiacal light as a function 
of elongation e can be approximated sufficiently by a power function of the 
form 

-K 

€ 
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RADIANCE ( 10 MAG. STARS/O 0 ) 



ELONGATION 


Figure 27. Collection of representative results between 25 and 110 deg 
for the radiance of the zodiacal light in the ecliptic [40, 27] . 




for e < .80 deg. Table 9 lists the values of the gradient exponent k for the- 
data of several observers. 

Table 9 shows that considerable discrepancies exist in the relative 
variation of the radiance along the ecliptic ( as obtained from different observa- 
tions) . The influence of the atmosphere on these data has yet to be determined. 

%■ 

TABLE 9. VALUES OF THE RADIANCE GRADIENT EXPONENT k [42, 43] 


Observer 

Range 

(deg) 

k Value 
(A., nm) 

Divari/Krylova [44] 

\ 

35 to 65 

2.4 (520.0) 

3.5 (520,0) 

Blackwell/Ingham [ 9] 

20 to 70 

2.4 (620.0) 

E Ivey/Roach [5] 

40 to 70 

1.9 (450.0) 

Robley [45] 

35 to 65 

2.3 (463.0) 
2.3 (616.5) 

Weinberg [46, 47, 48]" 

30 to 60 

2.1 (530.0), 

Peterson [49] 

25 to 60 

2.2 (435.5) 
2.2 (542.5) 
2.2 (638.0) 

Behr/Siedentopf [8] 

35 to 60 

1.5 (543.0) 

Divari/Krylova [50 ] 

40 to 80 

2.5 (540.0) 

Elsasser [51] 1 

30 to 60 

1.5 (visual) 

Fesenkov [34] 

. 30 to 60 

2.9 (-) 

Gillett [27] 

30 to 60 

2.1 (350.0- 
500.0) 

Smith/Roach/Owen [19] 

30 to 60 

2.2 (530.0) 

Regener/Van de Noord [52] 

25 to 35 

3.0 (435.0) 
3.0 (550.0) 


Figure 28 presents the distribution of the degree of polarization in the 
plane of the ecliptic [ 40] . Considerable discrepancies also exist in the values 
of the degree of polarization. For an elongation of 60 deg, the values shown 
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Figure 28. Collection of representative results between 20 and 160 deg 
for the polarization of the zodiacal light [40, 27] . 




vary from 0. 15 to 0.33. Generally, however, the polarization seems to 
increase when the elongation increases from 30 to ~ 65 deg. In the 50 to 
70-deg region, the degree of polarization changes little about the inflexion 
point. At elongations larger than 60 deg, the degree of. polarization decreases 
with increasing elongation. The polarization vector in the ecliptic, is dominated 
by the vertical component; i.e. , the component perpendicular to the ecliptic. 

i ' » - . 

In Table 10, the radiance and polarization values obtained by 
Weinberg [47] in the ecliptic are given for reference since his data were 
chosen as the base which the models are, compared. • The zodiacal light was’ 
observed at 530 nm with a photoelectric polarimeter atop Mt. Haleakala, 

Hawaii (geographic latitude 20° 43’ N, altitude 3052 m) between November 
1961 and May 1962. Simultaneous observations of the 557.7 nm airglow 
emission line were used toy-separate the components of night-ysky radiation. ; 
The field of view of polarimeter was 3. 2 deg in diameter. ' , 

• .. ' ’ ■ ’l J 

For radiance and polarization values within 30 deg of the sun, the data 
of Blackwell and Ingham [ 53] are given in Table 11. , i "■ 

The results of Smith, Roach, the Owen [19] for the variation of the 
radiance out of the ecliptic are given in Table 12. The data were obtained 
from Mt. Haleakala using a photoelectric photometer having a 4-deg field 
of view and 530 nm filter. The results are also presented in Figure 29 in an 
isophotal map where the circumference represents the plane of the' ecliptic 
with values of differential longitude, (X - increasing from 0 to 180 deg. 

The ecliptic latitude, 9, increases from 0 deg at the ecliptic. to 90 deg in the 
center of the circle. The photometric perturbation at 65 deg from the sun at 
high ecliptic latitudes has not been confirmed by other observers [54, 55] . 

I Dumont's [56] data for the average radiance out of the ecliptic are. 
given in Figure 30 for comparison. The minimum in the radiance for these V: 
data is at 70 deg and not at the ecliptic pole (90 deg) as indicated by Smith, 
et al. The difference between the data of Smith, et al. , and Dumont is a 
typical example of the discrepancy between different observers. The dis- 
crepancies between observers are again shown when a comparison of the 
radiance value of the zodiacal light at the ecliptic pole is made. Table 13 
presents data for such a comparison. Few observers have reported the 
polarization values at the ecliptic poles. Beggs et al. [57, 58] and Weinberg 
[54] found the values of ^ 20 and = 10 percent, respectively (Table 13). 

Color measurements have yielded color indices'both redder and bluer than 
sunlight (Table 14) . Also, the color index has been found to increase. [ 8] , 
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TABLE 12. ZODIACAL LIGHT IN ECLIPTIC COORDINATES [ 19] 


\ ■ 

■a 

■ 

Efl 

15 

20 

m 

30 

35 

40 

45 

50 

55 

60 

65 

WM 

HI 

80 

85 

180 

205 

195 

185 

175 

165 

155 

150 

140 

125 

130 

125 

125 

125 

125 

115 

115 

115 


175 

195 

190 

175 

165 

160 

150 

150 

135 

125 

130 

125 

125 

125 

125 

115 

115 

115 


170 

185 

175 

170 

160 

155 

150 

145 

135 

125 

130 

120 

125 

125 

125 

115 

115 

115 


165 

175 

170 

165 

160 

150 

145 

145 

135 

125 

130 

120 

125 

125 

125 

115 

115 

115 


160 

175 

170 

. 160 

155 

145 

145 

145 

135 

125 

130 

125 

120 

125 

125 

115 

115 

115 


155 

170 

165 

160 

155 

150 

145 

145 

135 

125 

130 

125 

125 

125 

125 

115 

115 

115 


150 

170 

165 

160 

155 

150 

150 

145 

140 

125 

130 

125 

125 

125 

125 

115 

115 

110 


145 

175 

170 

165 

155 

155 

155 

145 

140 

125 

135 

130 

125 

125 

125 

115 

115 

110 


140 

175 

170 

165 

160 

155 

155 

145 

145 

130 

135 

130 

120 

125 

125 

115 

115 

115 


135 

175 

170 

170 

160 

160 

155 

145 

145 

135 

135 

135 

120 

125 

125 

115 

115 

115 


130 

175 

175 

170 

165 

160 

155 

145 

145 

135 

135 

135 

.125 

125 

125 

115 

115 

115 


125 

180 

175 

175 

165 

165 

160 

145 

145 

135 

135 

135 

125 

125 

125 

115 

115 

115 


120 

185 

180 

180 

170 

165 

160 

145 

150 

135 

135 

135 

125 

125 

125 

120 

115 

115 


115 

190 

185 

180 

175 

170 

165 

150 

150 

135 

135 

135 

125 

125 

120 

,120 

115 

115 


110 

200 

190 

185 

175 

170 

165 

. 150 

150 

140 

135 

135 

125 

125 

120 

.120 

115 

115 


105 

210 

200 

190 

180 

. 175 

170 

150 

150 

145 

135 

135 

130 

125 

120 

120 

115 

115 


100 

220 

205 

195 

185 

175 

170 

155 

150 

145 

135 

135 

130 

125 

120 

120 

115 

115 


95 

235 

215 

205 

190 

180 

170 

155 

155 

150 

135 

135 

130 

125 

125 

120 

120 

120 


90 

250 

235 

215 

200 

185 

175 

155 

160 

150 

135 

135 

130 

125 

125 

125 

125 

125 


85 

265 

255 

235 

215 

195 

180 

155 

165 

155 

135 

140 

130 

125 

125 

125 

125 

130 


80 

290 

275 

255 

230 

210 

190 

165 

170 

155 

135 

140 

130 

125 

125 

125 

125 

130 

110 

75 

320 

300 

285 

245 

225 

200 

185 

175 

160 

135 

145 

130 

125 

130 

125 

130 

130 

110 

70 

(370) 

335 

310 

265 

240 

210 

205 

185 

175 

145 

145 

130 

135 

135 

130 

135 

130 

115 

65 

t 

(425) 

^70 

345 

285 

255 

225 

205 

205 

200 

170 

150 

140 

140 

140 

130 

135 

130 

115 

60 

(500) 

425 

380 

. 305 

270 

240 

215 

205 

205 

195 

165 

145 

145 

145 

135 

135 

130 

115 

55 

(600) 

505 

430 

330 

290 

255 

225 

205 

205 

200 

180 

160 

155 

150 

135 

135 

135 

120 

50 

(745) 

610 

485 

370 

315 

285 

230- 

205 

205 

195 

205 

170 

160 

150 

140 

140 

135 

120 

45 

(945) 

735 

555 

425 

345 

315 

240 

215 

210 

190 

205 

190 

170 

155 

140 

140 

135 

120 

40 

(1185) 

870 

655 

505 

390 

360 

250 

220 

215 

190 

195 

205 

170 

160 

140 

140 

135 

125 

35 

(1595) 

(1145) 

845 

690 

450 

410 

270 

225 

215 

190 

185 

200 

175 

165 

145 

145 

140 

. 125 

30 

(2330) 



905 

545 

465 

300 

230 

220 

195 

180 

195 

185 

170 

145 

150 

145 

125 

25 






520 


245 

225 

200 

180 

190 

190 

175 

150 
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145 

130 

20 









230 
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180 

185 
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180 

' 150 

150 

145 
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180 
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185 
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150 
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180 
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185 
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155 
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175 


185 

155 

160 



0 




(6150) 

(1820) 

(995) 

(525) 

(330) 

(250) 

(215) 

(180) 

(175) 

(185) 

(185) 
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(150) 

(135) 

































TABLE 13. RADIANCE AND' POLARIZATION OF THE ZODIACAL 
. • T . LIGHT AT THE ECLIPTIC POLE [59, 60] 


Observer 

(1) 

Pole- 

[Sio (vis) ] 

■ e = 4^ on the 
Ecliptic 
[Sio (vis)] 

Ratio ;~|j 

Polariza- 
tion at 
Pole 

'Elvey/Roach [ 5] 

73 

>'■ 1323 ... 

; 18.4 -A 


Beggs/etal £. [ 57 , 58] 

177 

66i y'. 

;3.73 ; ^ 

• . '•*' : * ‘ -* 

.. °- 20 

Weinberg [ 54] - 

186 j ‘ 

' 897 

■ ' 4. 82'' • •' 

0.09 

Smith/Roach/Owen [ 19] 

i i.0 Vi 

J 945 

8.60 ’ 


Dumont. [61] 

70 

i 

700 • 

- 10.0 

• *:» . > 

0.18 

Wolstencroft/Rose [62] 

' X143 -.; - 

840 ; 

5.87 

0.19 

Gillett [ 63] 
Sparrow/Ney [ 60] 

" 48.- 

50 . 

• ’■ v i i ■ ' 

. 416 

8.67 

0.21 

0.21 

• • ( 



TABLE 14. VALUES OF THE COLOR INDEX (Cl) OF THE 
ZODIACAL LIGHT [42] 


; — ! 

Observer 

Cl (zl) 

. Cl ( sun) 

<1 

x. A 

b g 

Behr/Siedentopf [8] 

0.63 . 

0.43 

+0.20 

444/542 

Divari/Asaad [64] 

0.35 ' 

0. 45 

-0.10 

414/541 

Nikol'skii [cf 42] 

0.56 . 

. ■ 


414/540 

Peterson [49] 

0.48 

' 0.45 

■ +0. 03 

435/543 

Divari/Krylora [ 50] 

’* 0.47 

0.63 

-0. 16 

406/543 

Robley [45] 

— 

— — 

+0.16 

463/528 

Divari/Krylova/Moroz [1121 

— 

.■ 

+0.14 

460/520 
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to decrease [64], and to remain constant [49]. If, indeed, an elongation 
dependence of the zodiacal light color index exists, it must be very small. 
However, if such a dependence is found, it will yield information on the 
dust-particle sizes. 

Figure 31 shows the relationship between the neutral-point position in 
elongation and the wavelength in the plane of the ecliptic [ 40] . For larger 
wavelengths, the neutral point tends to move closer to the sun. The data have 
a 25 -deg error spread so the linear nature shown may not be real. These 
data are important because neutral points calculated by models are dependent 
on the model chosen, at least more so than the general shape of the radiance 
curve produced by the models. 

The components of the zodiacal light caused by dust and electrons can 
be separated by the use of the spectrum of the zodiacal light. Blackwell and 
Ingham [9, 10, 11] determined that the electron density at 1 AU should be no 
greater than 120 electrons/cm 3 . The photoelectric observations of Beggs et al. 
[57, 58] which determine the electron contribution to be 16 ± 20 electrons/cm 3 
have been confirmed as to the order of magnitude from space vehicles [65] . 
Thus, it may be confidently stated that free electrons play a negligible role 
in the scattering of sunlight except in the corona. 


Photometric observations and recent results on the gegenschein have 
been summarized by Roosen [31, 66, 67] . He concludes that the gegenschein 
lies at the antisolar point, with no net displacement in either latitude or 
longitude to within 0.03 deg. The relative brightness of the gegenschein as a 
function of the distance in declination is shown in Figure 32. Its relative 
brightness as a function of longitude from the antisolar point is shown in 
Figure 33. The general shape of the gegenschein is oval. The shape param- 
eter, Q, is determined as a function of angular distance by the ratio of the 
brightness of the gegenschein and zodiacal light in a plane at right angles to 
the ecliptic to that in the plane of the ecliptic. Figure 34 presents the shape 
parameter as found by various observers. The discrepancy in the brightness 
of the gegenschein between various observers is given in Table 15. The ratio 

of the brightness of the gegenschein, B , to the brightness at the ecliptic 

G 


pole, B is given to find more consistency, but, again, a large amount of 

'scatter is evident since the maximum is 4.00, the minimum is 1.40, and the 
average is 2.10. The important result of Roosen is that to within 1 percent 


earth's shadow did not effect the brightness distribution of the gegenschein. 
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Figure 31. Relationship between the neutral -point position in elongation and the wavelength in the plane of the 
ecliptic. [The dashed line indicates the trend of Weinberg's data (error bars) for the dependence of the neu- 
tral point (zero polarization) in the ecliptic. The solid line represents result for the ice model P. 
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SMITH, ROACH & OWEN (1965) 5300A 
CHIPLONKAR & TILLU (1967) 4400A 
CHIPLONKAR & TILLU (1967) 5300A 
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TABLE 15. BRIGHTNESS OF THE GEGENSCHEIN, THE ECLIPTIC 
POLE, AND THEIR RATIO FOR VARIOUS OBSERVERS [ 66, 67] 



using V = 1.78 B (assuming solar color index) 



















SECTION VI. COMPARISON AND ANALYSIS OF 

MODEL CALCULATIONS ' ' 

A. Model Comparison 

From the preceding sections, one can see that to predict the observa- 
tional quantities of the zodiacal light, the* parameters of the models must be 
varied. Though the mathematical formulation for the problem is relatively 
simple, the observational uncertainties or lack of observational data has 
hindered the determination of definite values for the various parameters. In 
this section is presented an analytic approach to a method of variation of the 
parameters and a trial-and-error approach for matching the present data. 

A comparison of model predictions with observational data will give 
information on the range of the following: 

1. Particle spatial distribution. 

2. Particle size and numerical distribution. 

3. Particle shape and anisotropic effects. 

4. Particle composition (refractive index) . 

The following procedure is adopted in this section: . ,, 

1. Assuming i(0) is independent of Q and assuming that no radial 

dependence (m = 0) exists, it is shown that the radiance dependence on 

elongation e is approximately that which is observed. Also, assuming that 

the size distribution nice ) is a constant between a and a and other- 

min max 

wise zero, an indication of the number density parameter is given. As a 
further refinement to i(fl) , simple models using diffraction scattering theory 
are investigated. 

2. An indication of the material and a range on the particle size 

is found by comparing the radiance and polarization of the zodiacal light with 
that of single particles. 
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3. Further analysis in comparing polarization and radiance is 
accomplished by the use of scattering-intensity functions for a distribution 
over the particle size. This gives a further indication of particle size and 
materials. It is noted that a backscatter peak (gegenschein) in radiance is 
produced only by dielectric particles. 

, 4. Zodiacal light models are investigated for ecliptic variation of 

radiance and polarization. This follows the work by Giese [26]. First, a 
variation of materials is investigated, and then a variation of the parameters 
is investigated. Specific parameters considered are as follows: 

. . . . a. ; Spatial parameter (m) . ' - ‘ 

bi Number density variation parameter (p) . 

c. Number density constant (n^) . 

d. Minimum and maximum size parameter (a . and 

. min 

a ) . 

max 

e. . Material (refractive index m*) . 

f. Wavelength (\) . 

g. Ecliptic thickness parameter (R/W) . 

An analysis of the integration procedure and the influence of the maximum 
size cutoff is cursorily investigated. Finally, the various models are 
compared. 

5. Out-of-the-ecliptic calculations using the linear thickness model 
are presented; Thickness and distribution effects are investigated. 

6. Cylindrical and other shapes are considered as to their influence 
on the variations of the zodiacal light. 

Several conclusions have been made using restricted models for the 
zodiacal dust in the ecliptic. By comparing observations with models using 
spherical particles, it has been concluded that [54, 55] : 
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1. Only dielectric-component models produce a radiance enhance- 
ment corresponding to the gegenschein. 

2. Dielectric- component models produce an irregular decrease of 
the degree of polarization such as that observed at large elongations 160 

de g).- 

3. Dielectric-component models can give a maximum polarization 

at 70 deg. f - 

4. Metallic-component models cannot account for observed radiance 
and polarization. 

5. To include more than some tens of electrons per cm -3 at 1 AU, 
one must postulate the existence of a hypothetical nonpolarizing dust component. 

6. A simple power law for particle-size distribution and spatial 
distribution seem partially correct. 

7. The interplanetary matter has a low index of refraction 
(m =1.3- 1.8). 


B. Basic Heliocentric Model 


Results from the basic heliocentric model are presented to provide 
an understanding of the basic variation of radiance and polarization with 
particle number density and with elongation because of geometric considera- 
tions. The radiance equation for the heliocentric model is given by equation 
(33): 


r = 

X 


R E® 6\N 
X o 


. m+1 
sin e 



where 



_1 

N 


o 


a 

max 


/ 


a . 
min 


n 

o 




da 


j. . 
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and 


max 


N = f n / — \ 
o J o \ a / 
a . \ o / 

mm 


\ da 


An investigation was carried out in Section in. D with emphasis on the product 

1 

7T 

1 


I (e) = 
m . m+1 


J sin 0 d 0 


sin e e 


specifically, for m *= 0, 1, and 2. Employing 


Wo < £ > 


IT - £ 

sin e 


and 


I (e = 7r) = 1, 

m = 0 v ’ 


da 


one now uses that analysis assuming isotropic scattering (i. e. , = C = 

Q 

constant) and a distribution independent on the solar radial distant. Let 


i (0) + i (0) 2 max 


.o max / , \ 

xir / n o (an) 1:13 ’ 

o a . ' ' 


min 


then 


— \RE 5 A N 

L = *— 

X 47T 2 


o / l r-e 


ij + i 
1 o 


sin e 


(150) 


To use units of the number of 10th magnitude bolometric stars, one uses [69] 


solar constant = E 1. 39 x 10 3 w/m 2 


(m bol = 10 > 


star 


- — — I? 7 — 5 - = 0. 551 X 10 15 

2. 52 x 10 2 w/m 2 - 
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(151) 


where A = 500 nm. Using for — (ij + the value of 10 which is the average 

value for ice with an a on the order of 10 and 9 > 40 deg (i.e. , a ^ lju at 
500 nm) then at e = 40, ; 

I = ]L r~ L = 3. 8; ' 

m=0 sin e . , 

hence, to match Weinberg’ s data for e = 40 deg, the radiance equation (150) 
becomes ; c 

1.13 x 10 3 = 1.58 x 10 14 N ,(3. 8) (10) , ' 

. • - f ’ ; 

which gives , ' ^ ; . 

N ( e = 40) - i. 88 X io-» P arti f es i , ’ ■ : 

O V . ! : ' - cm' 3 ;, 

' r ' .i 

For e = 80 deg, one has , " . / ' i’ 

N (80)= 1.07xl0-» ‘ >artlC 3 lfeS * 

o v ' cm' 3 

and, for e =140 deg,' 

N (140) = 1.17 X 10-3 . 

o v ' cm 3 4 

The values for the radiance, using the various values of N q compared with 

Weinberg’s data are presented in Table 16 and Figure 35. The last column 
(A) is the difference between the third and sixth column which gives an 
indication of the average total intensity as a function of 9 , although the 
radial dependence is also, contained in this difference. 

The 30 to 60 deg exponential gradient for these models is K = 1.2, 
compared to 2. 1 from Weinberg’ s data. Hence, in the basic model, the 
radiance rises slower in the solar direction than measured data. The wave- 
length or spectral dependence is the same as the sun (i.e. , a GO-type star). 


to obtain the value of 


— = 158. 28X 10 -12 ■ 1 ° ma ?-. 3ters - cmi 


4t 


0 deg 
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Figure 35. Comparison of the basic heliocentric model 
with Weinberg' s data. 
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TABLE 16. COMPARISON OF THE BASIC HELIOCENTRIC MODEL 

WITH WEINBERG'S DATA 


e 

f 

o 

Weinberg 

s io( v ) 

Heliocentric Models 

A 

N Q (40) 

N q (80) 

N o (l40) 

30 

5.23 

2200 

1555 

886 

967 

1233 

40 

3.80 

1130 

1130 

644 

• 704 

426 

60 

2.41 

520 

717 

408 

446 

74 

80 

1.77 

300 

526 

300 

328 

-28 

100 

1.41 

215 

419 

239 

261 

-46 

120 

1.20 

185 

357 

203 

222 

-37 

140 

1.08 

200 

321 

183 

200 

0 

160 

1.02 

220 

303 

173 

189 

+31 


For a distribution of particle sizes of 

1 a . < a < a 

mm max 

0 otherwise , 

and for a . =0.5 /xm and a =1.5 |im, then 
mm max 

N = n (a . - a ) 
o o min max' 

or 

N 

_ o 
o l/i m 

Hence, for N q ( 140), the parameter n^ has the value 

n ( 140)= 1. 17 x 10~ 9 • 

o cm 



120 















dc . 

In conclusion, the radiance dependence on is not very great, as was noted 
in earlier theories [13] , although the lack of this dependence is quite noticeable. 


C. Basic Heliocentric Model with Spatial Dependence 


The radiance variation of the basic heliocentric model (HM) does not 
correspond to Weinberg’ s data. To improve the basic model, a spatial de- 
pendence for the number distribution of the form . 

s(r) a* l/r m 


is included in the model. From the analysis of I , it is evident that m - 1. 

m 

Hence, 



1 + cos e 
sin 2 e 


(152) 


will give e" 2 dependence for the radiance as e — 0. This is comparable with 
observed data (Table 9) . Although for I _ , we have a e “ 3 dependence as 

e — 0, e“ 2 dependence can be obtained for a higher range of e. The result of 

these radial dependence models, normalized to Weinberg' s value at e= 140 

deg, is shown in Table 17 and Figure 36. Comparison of the 30- to 60-deg 

exponential gradient gives k =1.7 and k = 2.6. 

m = 1 m — z 

TABLE 17. HELIOCENTRIC MODEL WITH SPATIAL DEPENDENCE 


e 

'm = 1 

*m = 2 

L; ^ 

m = 1 

L b 

m = 2 

30 

7.46 

12.2 

2635 

6310 

40 

4.27 

5. 52 

1508 

2855 

60 

2.00 

1.94 

. 706 

1050 

80 

1.21 

1.00 

427 

517 

100 

0.852 

0.641 

301 

332 

120 

0.666 

0.4(80) 

235 

248 

140 

0.566. 

0.387 

200 

200 

160 

0.515 

.0.345 

182 

178 


a. N = 3. 53 x 10“ 13 

o 

b. N = 5. 17 x 10 -13 

o 


M454 
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D. Belt Model for Isotropic Scattering 


The variation of radiance with respect to elongation can be improved 
by the use of the belt model which is described in Section III. For simplicity, 
assume an isotropic scattering intensity functions; i. e. , I 


Me) +i (e) 


= constant 


and m = 0. As before, the constant is chosen to be 10. In units of 10th 
magnitude stars per square degree, the radiance equation (59) becomes 




158.28 x IQ 12 X 10 
sin e 


Z 


j = i 


V 1 

N. / de 


0 


J 


where 


0 . = sin/— sin e \ 

J V) J 


The values (R/r.) 


-l 


are chosen to be 0. 723 and 1.00; i.e. , the orbits of 


Venus and Earth. For sufficiently small elongations, the line of sight crosses 

the Venus-Earth belt twice. Note that N. and/or 0 . could be varied, but be- 

• / 3 3 

cause of the dependence of 0^.-' on 'the elongation, one determines, for ease of 

3 • 

calculation, the 0.’s by choosing- the r ’s. 

j . ; j 

' For :NZ ? in units of 10" 1 , 3 particles/cm 3 (using Weinberg’ s data), the 
following conditions exist: For e = 40 deg, 


1130 = o^279 [N 3 (l. 09528 - 0.69813) 

+ N 2 (2.04631 - 1.09528) 

+ N 3 (2.44349 - 2.04631) 

' + N 4 (3.14159 - 2, 44349)1. 
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For e = 60 deg: 


158.28 
0. 86603 


[N 3 (2. 09439 - 1.04720) + N 4 (3. 14159 - 2.09439)] . 


For e = 120 deg: 

'“' i ^'^ 3 - 14159 - 2 - 09440 ) 1 - 


N 2 represents the density inside r = 0. 703, N 3 represents the density between 
r = 0. 703 and 1. 00, and N 4 represents the density for r > 1.00. These equa- 
tions yield 


N 2 = 2.653 x 10" 13 
N 3 = 1.747 x 10 -13 
N 4 = 0.969 x 10“ 13 


particles/em 3 


where the factor 10“ 13 has been included. The comparison of the belt model 
and Weinberg’ s data is given in Table 18. Also, the comparison is made in 
Figure 37 which shows the basic heliocentric model for m = 0. The results 
are promising although more parameters are introduced relatively ad hoc. 
The 30- to 60 -deg exponential gradient for the belt model is k = 1.8. 


The density calculated decreases with distance for the three separate 
concentric regions. Since i(0) is constant, the effect could represent either 
the actual spatial dependence or the diffraction peak, or both. The diffraction 
phenomena will be investigated in Section F. 


In the belt model, when r. = 1. 524 and e = 170 deg, then only approxi- 


mately the last 10 deg of the scattering intensity functions for particles outside 
Mars orbit. contribute. Thus, particles with very large backscatter could be 
the main contributors to gegenschein. This effect could be a type of opposi- 
tion effect for large particles (i.e. , with a radius greater than*! cm) outside 
the Martian orbit. 


E. Approximation 

In the calculation of scattered light, approximations can be used to 
represent the angular distribution of the scattered light. These approximations 
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Figure 37. Comparison of a belt model with Weinberg' s data. 
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TABLE 18. COMPARISON OF A BELT MODEL WITH WEINBERG' S DATA 


€ 

Weinberg 

A Belt Model 

30 

2200 

1782 

f 40 

1130 , 

1130 

50 

730 

678 

60 

520 

520 

70 

385 

404 { 

r. 

.80 

300 

315 j .. 

90 

250 

240 

100 

215 

217 : “ 

110 

195 

199 

120 

185 

185 

130 

190 

174 

140 

200 

166 

150 

205 

160 

160 

220 

156 : 


were empirical and usually fit a particular scattering curve. The advantage of 
the approximations over exact expressions is easily seen when the labor 
needed to compute the exact scattering functions is corhpared to that needed for 
an approximation. Some of the earlier approximations to scattering functions 
of particles were those used to explain diffusely scattered. light from a planet. 
As may be expected, no diffraction peaks were included. Figure 38 shows four 
such functions [2] : 
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a. 


sin 0-0 cos 0 


b. 


c. 


d. 


7T 


(Lambert) 


. /0 - 7T\ /0 _ 7T \ . /0 - 7r\ 

1 - sin v~ 2 / 13,1 Hr ~ ) 1 n cot Hr - / 

( Lommel-Seeliger) 

/ff _ a \ 2 (VanRhijn) 

\ ir / 150 ^ 0 ^ 180 deg 

0.655 - 0.711 sin 0 - 0.345 cos 0 + 0.406 sin 3 0 


(Van Rhijn) 
40 deg ^ 0 


Van Rhijn found his phase function from a best fit to observed radiance curves 
for zodiacal light. Lambert’s law and the Lommel-Seeliger law have also 
been applied to such calculations [2] . 


To explain scattering from interplanetary media at small elongations 
from the sun, empirical scattering functions are needed which include a dif- 
fraction component and a reflection component. Allen [12] proposed that the 
F corona of the sun is due to diffraction on interplanetary material. He used 
an approximation to the diffraction function given as 


C a 2 

2 + a 3 sin 3 0 


(153) 


It is plotted in Figure 39, although it was not intended to be used out to 180 deg. 
At almost the same time. Van de Hulst [13] also explaining the F corona 
separated the scattering function into a diffracted and reflected part, using 
Lambert's law for the reflection term. This separation procedure is valid 
for large particles. Van de Hulst’ s function is 


J t 2 (o? sin 0 ) + co_ 
a 2 sin 2 0 47r 

where co is the albedo of the particles. Since the calculations of Van de Hulst 
represent a simple model of the zodiacal light which incorporates the diffrac- 
tion effects, it will be considered further in Section F. 

Schoenberg [2] used the general formula 


1 - p cos © + q. cos 2 0 


(154) 
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SCATTERING ANGLE (DEG) 

Figure 39a. Approximate scattering functions with a diffraction peak, 
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; 



SCATTERING ANGLE 

Figure 39b. The scattering functions of Gillett [ 27] . 
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to characterize the properties of scattering in the atmosphere. He used 
p = 2. 7 and q = 3. 0. With p = 0 and q = 1, this formula reduces to that needed 
for a Rayleigh scatterer. To correct for the effects of relatively large 
particles, V. A. Kratt modified the Schoenberg formula to [ 70] 


„ , -3 0 -37r/2> 2 « 

1 - p(e - e ) + q cos' 1 0 


(155) 


'This has been further modified for corona, etc. , by adding or subtracting 

terms [e.g. , 1 - p(e 3 9 -e 37r ^ 2 ) + q cos 2 0 + S e 36 9 ] [113]. Divari 
[71] used the formula 


-3 0 


17.9 


[ 1 + 22. 3 (e -0.009) + 0.02 cos 2 ©] 


(156) 


to calculate the contribution of a dust cloud around the earth to the radiance of 
zodiacal light. 

Banderman [36] , by splitting the scattering function into a diffracted 
and a reflected component, found that, according to his meteor models, the 
zodiacal light radiance is given by 

0.1 

2n 

The particle albedo is 0. 1 and the ratio of the mean particle radius to the mean 
square of the particle radius is about 0.264. The wavelength used for the cal- 
culations was 0. 53 p. 

<► 

• .Gillett [27] , also attempting to model the zodiacal light for different 
models of the. cloud of interplanetary medium, used as a best fit 

| 3/4 cos 2 © 0 ^ 7r/2 

— -- - sin 0 + (158) 

sin — ( 0 0 > n/2 

in his constant thickness models (Paragraph E of Section Hi) . 

\ 

Each function discussed above is plotted in Figure 39 and has been 
normalized to one at 0 = 90 deg. Apparently, Bandermann* s phase function 


1 + 


(1.056V 


1 + e 


1.056tt 


1.056 0 1 , , 

e +- (0.14) (0 7r) 


-3 


(157) 
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is the best approximation to the shape of the radiance curve for zodiacal light. 
Table 19 lists the functions and their authors. 

TABLE 19. PHASE FUNCTION, APPROXIMATIONS 



Author 

* 

F. Schoenberg [2] 

V. A. Kratt [70] 

■ • _;o' 

N. Bi Divari [71] 

C. W. Allen [12] 

H. C. Van de Hulst [ 13] 

L. W. Banderman [36] 
F. C. Gillett [27] 

Van Rhijn [2] 

£. 

Lommel-Seeliger [2] 
Lambert [2] 


Function 


1 + p cos 9 + q cos 2 0 

1 / —3© — 3 j 2 •. o r\ 

1 + p(e -e ) + q cos' 1 0 


1 30 

— —[1 + 22.3 (e -0.009)+ 0.02 cos 2 0] 

1 I • >7 

C a 2 /(2 +'oi 3 sin 3 0) 

Jj 2 (a sin0)/o; 2 sin 2 0 + co/4tt 

T 1 e 1 - O56 0 + O.O7 (0*)-* 

2ir , 1.056tt ' ’ 

1 + e 

^ • 3/4 cos 2 0 0^ 7r/2 

— - +-sin:0 + - 

sin-?- 0 0 — 7r/2 


0.655 - 0.711 sin 0 - 0.345 cos 0 + 0.406 sin°0 

(~ ) 2 

1 - sln (^T 2- ) tan (^T 1 ) ta cot (■ § T 1 ) 

sin 0 / 7r - 9 cos 0 /tt 


Obviously, many approximate scattering functions exist which can be 
used. For zodiacal light, the most common are given as a sum of a reflected 
and a diffracted component. Chandra shekhar [114] points out that any phase 
function can be represented by a series of Legendre polynomials 
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\ 


h(e) = £ u n p n (cos 0) . 

n 

co is a constant for a given particle size and material. In general, u^isa 

function of a and m. Hunt [72] uses a discrete least-squares method to find- 
the phase function found from Mie theory calculations as a series of Legendre 
polynomials. Such representations of the phase function were found for both 
single scattering and scattering from a polydispersion of spherical particles, 
specifically the haze and cloud models of Diermendjian [115] were used. 
Hunt's method gives an overall relative error of 10 -5 for single particles and 
1CT 8 for the distributions. It should be noted that phase functions for po ^dis- 
persed collections of particles are independent of the density of particles due 
to the. normalization constant. 

F. Diffraction - Reflection Models [12, 13] 

Van de Hulst [ 37] , in a study of zodiacal light in the solar corona, 
separated the differential cross section into a diffraction component and a 
reflection component 


do _ / do \ / da \ 

d« " ^ dSI ) , + \ dS2 ) 

d r 


(159) 


which is a first-order approximation that is valid for particles with a size 
parameter greater than 5. The first term is the differential scattering 
caused by forward diffraction of a sphere. It is dependent of the particle 
surface structure and refractive index. ( Fraunhofer diffraction gives 



J (ka sin 0 ) 
ka sin 0 



j2 (k a 0 ) 



( 160 ) 


The second term is for reflection of a diffusely reflecting sphere 
(Lambertian), 




(sin 0 - © cos 0 ) , 


( 161 ) 
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which includes an albedo, factor y [ 13] . The reflection differential cross 
section is independent of the incident wavelength. " These scattering functions 
are shown in Figure 40 for a = 10 and y = 0. 33 f 73] . . . I 

From the photometric observation of Richter [74, 75] , scattering from 
meteorites of 10 2 to 10 _1 cm in diameter approximately agrees with Van de 
Hulst’s equation (159). In the range of 10 2 cm, little variation of radiance 
exists with the scattering angle outside the forward diffraction zone. This 
gives rise to a simpler form for the reflection differential cross section. 
Assume the total cross section is the geometric cross section: 


= / Cl) dfl = f — dfl = 7ra 2 , 


k 2 


then for an isotropic scattering, i.e., ij = = constant, 


2 (*1 4 * 2 ) 

w 


4tt = Tra 2 ; 


hence, for a sphere of radius at, 

•5 ‘ . 

h ^ 

(W ZkF { 1 z) 4 : * 


An albedo factor y is added to give 

• . * , i 

,2 


/ da \ _ ya 2 
Vdfl/ ~ 4 * 


■i 


(162) 


for which Figure 41 displays the total differential cross section. The approxi- 
mation is also good forj metallic spheres [76] . Following Van de Hulst [13] 
and using equation (33) , the radiance for fixed wavelength and particle size 
is given by : ;; 


R 6A N >a 2 ' 


L x - 


4e ' 


’ F(e) 


where 


F(e) = 


sin e 


7T 

5 4tt 


1 do - (9,) 


d0 = - 


4e 


a 2 sin e 


/(&) 


d0 


(163) 



DIFFERENTIAL CROSS SECTION (cr/) 



Figure 40. The differential cross section for a Lambertian sphere, the 
forward diffraction peak for a sphere, and the sum 
(a = 10, A = 0. 5 pm, y *= 0.33) . 
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DIFFERENTIAL CROSS SECTION (. 


N 

E 


SCATTERING ANGLE 

Figure 41. The differential cross section for an isotropic .scatterer, the 
the forward diffraction peak for a sphere, and their sum 
(o' = 10, \ = 0. 5 nm, y* 0. 33) . 




T“ ' • 

and where the particle distribution has an r -1 radial dependence. For isotropic 
scatter ers with y = 1, the reflection part becomes 


F 

r 


(€)■ = 


eU - e) 

sin e 


(164) 


which is approximately constant, varying from tt 2 /4 at e = 7r/2 to 0 at e = 0 
or 7T. For the Lambertian function, 


F (e) = y — — — (2 + 2 cos e + e sin e) 
r w . sine 3f ' ' 


(165) 


Both forms of the reflection functions remain relatively constant, giving 
radiance an e -1 dependence. 

For the diffraction component, Van de Hulst obtained 

F^(e) = ka ip (ka e) , (166) 


2 

( 2 .. Ju L z .l ) dz . (1 67 ) 

Since the diffraction term is only accurate for small values of 0 , the factor 
e/sin e has been omitted and the range of integration has been extended to °o . 
The values of ip(z) are given in Table 20. 

For large values of z, ip(z) fluctuates around 2/ttz 2 [ 77] . The radiance 
diffracted by large particles is therefore proportional to e -3 . The values of 
F d (e) calculated by Van de Hulst [ 13] are shown in Figure 42. The radiance 
variation can be divided into ttiree regions. 

Region Radiance Dependence 


(Lambertian) 



where 

00 

H z ) = I 

z 
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I. 


5°.7 10 ° 30 * 60 ° 90 ° 180 ° 


Figure 42. The scattering dependence F d (e) for diffraction of spherical 

particles for a = 10, 100, and 1000 {curves. (a) and (b) are for a 
isotropic scatterer and a Lambertian sphere (y = 0.1), respectively [ 13] } . 

TABLE 20. ip( z) VALUES 


1.698= 16/3 tt 


0.0156 


0.779 


0.0127. 


0.225 


0.0112 


0. 0074 


0,040 

0.030 


0.0061 


5 


t 















Comparison with' Weinberg' s data 5 is given in Table 2T and Figure 43 for N q 
established at e= 140 deg and for q; =T } 10, and 100. Although using 
a;=l for the diffraction peak is unrealistic, it gives the. best fit and also 
gives an mdication of’the form of da/dfi . 

As early as 1893, Searle, in modeling the zodiacal light by a dust cloud, 
used various simple scattering functions and a density distribution proportional 
to a constant or l/r. It is interesting to compare his and Van de Hulst' s 
'results, although Searle* s models do not include the diffraction contribution 
which originated with Van de Hulst* s work. 

Searle’ s models are: ’v 

Model Scatterihg Function Form Spatial Distribution Form 

■ I , Sin 0-9 cos 0 Constant 

i,..) .- ]: 'f < *■••••• (Lambertian) ' ' ■ 


, : n. : j 1 - cos 0 4 » Constant . . , 

lv IE 1- cos 0 1 l/r 

IV 9 Constant 

b. *X\\ ' ■ V 1 . j 

V 0 l/r 

■ VI 1 0 2 \ J l’/r '' 

The scattering distributions reflect observations made of the asteroid scatter- 
ing functions. Searle’ s [ l] results of the normalized values are shown in 
Table 22. It is worthwhile to quote this early work:' ‘ 

:"A remark obviously suggested by the inspection of Table XXVII .[Table 
;22 in this report] is that for small values of e the hypothesis relating to 
the density of the meteoric 1 matter is more important than that relating 
to the phases of the separate particles, while for large values of e this 
Hast hypothesis chiefly determines the result.' Compare, for example, 
the values derived frdm Formulas HI and V,' both of which assume the 
density to vary inversely-to r. For small values of ey these’ values are 
practically the same, although in the first case the projection of the 
illuminated surface, and in the second that surface itself, is assumed 
to represent the quantity of light reflected. Even when we employ the 
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TABLE 21. DIFFRACTION MODEL CALCULATIONS 
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TABLE 22. COMPARISON OF HYPOTHESES [ l] 
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very forced hypothesis of Formula VI respecting the effect of small 
phases but retain the same hypothesis respecting density, we reduce 
the light only by half a magnitude. On the other hand, for the same 
small values of e, Formulas II and IV, in which the hypothesis of 
uniform density is adopted, agree with each other and disagree with 
the remaining formulas. But for large values of e, according to the 
respective hypotheses adopted for the phases of the particles, the first 
three formulas agree with each other and disagree with Formulas IV, 
V, and VI, in which the minimum of light occurs before opposition is 
reached, while the light at opposition becomes equal to that at an 
elongation of agout 85 deg to about 125 deg." 

"All the formulas agree in representing fairly well the rapid 
diminution in brightness of the zodiacal light at small elongations 
and its much slighter subsequent changes, while the fourth and fifth 
. formulas show that a representation of the phenomenon of gegenschein. 

can be theoretically obtained even when the hypothesis respecting the 
1 phases requires very moderate assumptions, as will appear more 
clearly below. Upon the whole, the fifth formula best expresses the 
observed facts; the sixth unduly magnifies the relative amount of light 
' at opposition, although so faint a' gegenschein as that indicated by the ■ 
fifth formula would possibly never be discovered." 

This early work describes the important features of a model and shows 
how relatively simple models can depict the dominant effects on the zodiacal , 
radiance. 

Roosen [31, 66, 67] employs a reflection model for the scattering 
functions to predict the radiance of the gegenschein. Roosen was able to fit 
Weinberg* s observations of the gegenschein by using a 0. 023 mag/ deg phase 
function and a spatial distribution corresponding to that for asteroids (Fig. 
44a) . He concludes from the observed lack of effect from the earth* s shadow 
on illumination of the gegenschein that a spatial distribution of particles which 
follows an inverse power law cannot satisfactorily explain the gegenschein 
(Figs; 44b and 44c) and that there must be an increase in the spatial density 
outside the earth* s orbit. He suggests the asteroid belt is the source for 
this material. The spatial distribution of asteroids gives support to this 
argument. Recent observations at 500 nm from OSO-6 [116] show, based 
on a preliminary calibration, a gradient for the gegenschein radiance which 
is in accord with Roosen* s data (Fig. 44d) . 
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Figure 44a. The asteroidal phase function determined by Gehrels [ 78] 
showing the opposition effect. The phase function with the 
opposition effect removed [66, 67, 79] . 


G. Polarization - Empirical 

The basic heliocentric models have been used with a simple scattering 
intensity function and therefore exhibited no polarization. In this section we 
consider the experimental investigations of Richter [14, 74, 75] on the polar- 
ization from single particle scattering and distributions of particles. An 
investigation on the polarization attributed to several particle sizes and 
materials will give an indication of the material and particle- size distribution 
of the interplanetary dust. 
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Figure 44c. Gegenschein brightnesses for models assuming the material is 

distributed according to an inverse power law R m (where R is the helio- 
centric distance) for an isotropic phase function [66, 67, 79] . 
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To introduce polarization in simple zodiacal light models, Ingham [53] 
used diffraction scattering and the empirical form 


p(e) 


sin 2 © 

q - sin 2 (0 ) 


(168) 


for the degree of polarization, q is a constant depending on the nature of the 
material. Figure 45 plots p(0 ) versus the scattering angle 0 for q= 3, 4, 5, 
6, and 11. Note that p(0 ) given by equation (168) is symmetric around 
0 = 90 deg and that the maximum value of polarization is given by 


maxp(0) =p(0= 90 deg) = 

An inversion gives: 

, . i 
^ p ( © = 90 deg) 

The justification for assuming a form such as this results from the labora- ? 
tory studies of Richter [74, 75, 80] and Kloverstrom and Rense [73] on 
photometric properties of meteorites. 

Figures 46 and 47 show a portion of Richter' s results. In Figure 46a, 
radiance dependence on scattering angle is shown for pulverized iron, graphite, 
and quartz dust in the =* 1- to 10-Mm size range. The size distributions were 
roughly Gaussian. The modes are indicated in Figures 47a, b, and c which 
show the variation of polarization with angle. According to Little [86] , 
Richter's results show that dielectric and metallic powders differ from one 
another in the observed radiance (directly related to the scattering intensity 
function) . The radiance of the dielectrics decreased mbre or less smoothly 
to 0 = 140 deg and then rose again slightly toward 180 deg; the metallic 
variation declined to about 100 deg, rose to a secondary maximum near 150 
deg, and then declined toward 180 deg. The polarization curves of iron and 
quartz powder are similar in that both rise to a flat maximum not exceeding 
20 percent near 0 = 90 deg (Fig. 47), giving experimental foundation for 
Ingham's polarization function. Quartz powder polarization differs from the 
iron powder since the maximum is shifted towards larger values of 0 (Fig. 

47b) . Graphite dust shows a difference in that the observed polarization 
attains a maximum of 30 percent near 55 deg. 
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Figure 45. The degree of polarization' predicted by Ingham’ s empirical form for 

q = 3, 4, 5, 6, and 11. 


RADIANCE (ARB. UNITS) 



Figure 46a. Radiance of powders with micron-size particles for 
iron, graphite, and quartz [80] . =•, 



Figure 46b. Radiance of 6-cm iron meteorite compared with the 


Lambertian curve [75] . 











Mie theory (spherical particles) can be used to match the radiance and 
polarization of the irregular bodies of the powder studied by Richter. The 
results [86] indicate that dielectrics with a small imaginary component in the 
refractive index in the Mie theory can represent Richter’s data, although the 
radiance and polarization of the powder in Richter’s data were not reproduced 
by Mie theory. 

In Figures 46b and 47b, the observed radiance and polarization for a 
6-cm iron meteorite is given for comparison with that from the micron 
powders. The radiance does not show the diffraction peak but only the reflec- 
tion component. The 47-percent polarization maximum occurs near 40 deg. 

In general, the stony meteorites have a smaller radiance and polarization than 
metallic meteorites [74, 75]. For submicron particles, Rayleigh’s theory 
can be applied. 

Following Ingham [53] , the heliocentric model is used with an isotropic 
scattering function: 



and an empirical polarization function 
i - i 

p(0 ) = ~r r- = sin 2 0 /(q - sin 2 0 ) 

' ' i, + i 

l 2 

Hence, 
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where 


a 


max 


g (9) = / [ii(9» a ) - h (Q* <*■)! 


n(a ) , 

— * — L d a 


k2 


a 


min 


V £) - 


7T =-j— 

r . ni da 

/ SU1 0 ‘ST de 


sin e e 
a 


da 

dfi 


r 

j 

a 


max 


(ij (©, a ) + i 2 (©.a )) da 


mm 


Now, assuming n(a ) = a P , the polarization of the zodiacal light 'is expressed 
by 


7r . 2+m 
sin 


/ 


q - sim 0 


d0 


p x <£ > = ir 


f sin m 0 d0 


(169) 


with 


V £ > 


yC \ 2 


, „ o . m+1 

16 it sin e 


a 

J d0 sin m 0 J 


max 

a P da . (170) 


a 


min 


Ingham* s results for m = 1 q = 3. 2 and m = 3/2 q = 3. 3 for p(0 ) and p(e) are'\ 
given in Table 23. Figure 48 gives the degree of polarization versus elongation 
for q= 3, 4, and 5 for three values of m(0, 1, 2). Ingham's choice of q was 
arrived at by maintaining the electron density derived by Blackwell and Ingham 
[9, 10, 11] and Giese [26] . The q*s give values of P(90 deg) which are higher 
than any P(0 ) measured by Richter for freely suspended particles. Dollfus 
[ 81] , however, measures a maximum polarization greater than this at angles 
greater than 90 deg. However, the powders of Dollfus were spread on a sur- 
face and were not freely suspended. 
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TABLE 23. INGHAM’S POLARIZATION RESULTS FOR THE 
ISOTROPIC SCATTERS [53] 


e 

P(0) = 

■ 

sin 2 0 



P(e) 

P(e) 

q - sin 2 0 


q= 3.2 

q = 3. 3 

m = 1', q = 3. 2 

m = 3/2, q = 3.3 

0 

0.000 

0.000 

0.279 

0.290 

10 

0.009 

0.009 

0.281 

0.290 

20 

0.038 

0.037 

0.287 

0.295 

30 

0.085 

0.082 

0.296 

0.300 

40 

0.148 

0.143 

0.306 

0.309 

50 

0.225 

0.216 

0.315 

0.317 

1 60 

0.306 

0.294 

0.320 

0. 324 

70 

0.381 

0.365 

0.317 

0.324 

80 

0.435 

0.416 

0.303 

0.312 

90 

0.455 

0.435 

0.279 

0.290 

100 

0.435 

0.416 

0.246 

0.255 

110 

0.381 

0.365 

0.203 

0.213 

120 

0.306 

0.294 

0.159 

0.169 

130 

0.225 

0.216 

( 0.115 

0.123 

140 

0.148 

0.143 

0. 076 

0.082 

150 

0.085 

0.082 

0.045 

0.051 

160 

0.038 

0.037 

' 0. 018 

0.029 

170 

0.009 

0.009 

0.007 

0. 003 

180 

0.000 

- 0.000 

0.000 

0.000 
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: Figure 48. The degree of polarization for the zodiacal light for Ingham’s empirical model 

for q = 3, 4, and 5 and m = 0, 1, and 2. 
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H. Polarization-Single-Particle Scattering 

To reproduce the variation of radiance with elongation, the empirical 
scattering function used by Ingham does as well as more sophisticated ones, 
but polarization depends markedly upon the nature of the material and physical 
dimension of the particle. Mie theory can be used to investigate the polariza- 
tion and radiance dependence on particle size and material [24, 82, 83] . The 
results of such an investigation can be used to predict or isolate certain ranges 
of these variables for a given light- scattering situation. In this section the 
radiance and polarization from single particles is discussed. In Figures 49 
and 50 the total scattering intensity functions versus the scattering angle for a 
single spherical particle with a refraction index of 1. 33 (dielectric ice) is 
shown for various values of a (=27r a/X) . The corresponding polarizations are 
shown in Figures 51 and 52. For small particles, the scattering follows the 
Rayleigh law with a symmetric polarization of about 90 deg. As a increases, 
the diffraction peak appears, the patterns become asymmetrical, and the num- 
ber of maxima and minima in the scattering intensity functions increases. The 
polarization deviates from the Rayleigh case as a increases and oscillates about 
zero more with the higher a's. This can be seen in Figure 53 which shows the 
degree of polarization for a particle with a = 20. Figure 54 shows the total 
intensity for this particle. The diffraction peak has become very narrow and 
an increase in backward scattering occurs. 

Figures 55 and 56 show the dependence of the scattering intensity func- 
tion and polarization upon the absorption term of the refractive index. In these 
figures, m* = 1.33 + ik(k = 0.00, 0.05, 0.10, 0.50) and a = 5.0. The total 
scattering intensity function decreases by a factor of 2 in the forward direction, 
and 10 in the backward direction. The polarization for a small absorption term 
does not differ much from the pure dielectric, but for k = 0. 50 the polarization 
is totally positive, showing a metallic nature. The behavior is further shown 
in Figures 57, 58, 59, and 60, where data for dielectric and metallic particles 
of size parameter a = 5 are presented. The polarization from the dielectrics 
oscillates in sign while the metals give a predominantly positive polarization. 

An analysis using single-site particles indicates that the positive 
polarization of the zodiacal light can be obtained from small particles (a < l) 
or metallic particles with a greater than 1. This must be considered a 
primitive result because of the integration over particle size and scattering 
angle that is necessary in a model of zodiacal light. 
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Figure 50. Total scattering intensity functions for spherical ice particles 

(m* = 1.33) for a = 6 - 10. 
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Figure 53. The degree of polarization for a spherical ice particle 
- , (m* = 1. 33) for of = 20. 



Figure 54. The total intensity scattering function for a spherical ice 
particle (m* = 1.33) for a = 20. 
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Figure 55. Total scattering intensity function for spherical particles with 
various values of the absorptive part of the refractive, index (m* =1.33 - iK) . 
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Figure 57. Total scattering intensity functions for various spherical 
dielebtric particles with a = 5. 0. 







Figure 58. Degree of polarization for' the dielectric particles of Figure 57. 
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Figure 59. Total scattering intensity functions for various spherical 
metallic particles with a = 5. 0. 
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I. Scattering from Distributions 


A distribution of particles will tend to average and smooth the oscil- 
lations of the radiance and polarizations curves of single particles. The 
results using examples of each of -the four categories of distributions are 
compared and discussed in this section. Particle sizes between a = 1 and 
a = 20 and a = 0. 1 and a. = 20 were used, since this represents a possible 
size range of particles of zodiacal dust and does not require excessively long 
computer times. - 

i 


Power Law Distribution. Since this distribution is used extensively 
in zodiacal light calculations, it has been the most- frequently used in the . 
present studies of light scattering from a distribution of spherical particles, 
particularly with p = -4. The materials considered were water (ice) , quartz, 
titanite, carbon, iron, zinc, and nickel. The scattering intensity functions 
were calculated from 


i. 
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a 


max 

/ i (0, m* a) n( a) ■ do; 

a j 

min 
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max 

r 

J n(a 


OL 

min 


(a ) da 


j = 1, 2, 


(171) : i 

1 • ' 


where a . is the minimum particle size in the distribution and a is the* 
min max 

maximum particle size. i.(9 , m*, a) was calculated from the Mie theory. 

Figures 61a and 61b show the intensity versus scattering angle for the materials. 
(Each curve is displaced a cycle upward for clarity.) The backscatter peak 
is visible for the dielectrics, this being one of the well-known distinguishing 
characteristics between metals and dielectrics. Another distinguishing 
characteristic is the degree of polarization, being positive for metals in which 
the absorption term k is at least a fourth of the real part of the index of 
refraction [39] . Dielectrics produce negative polarizations at some scatter- 
ing angles. For the materials used, these properties are clearly visible in 
Figures 62 and 63. With ice, the radiance and polarization for distributions 

with a >20 does not vary much from that found with particles with a = 
max max 

20, as seen from Figures 64 and 65. For iron particles, the polarization is 

changed very little by particles with a > 5 as seen from Figures 66 and 67. 
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Figure 61a. Radiance versus scattering angle for dielectrics. 
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Figure 61b. Radiance versus scattering angle for metals. [Curves 
are displaced a cycle, n (a) ~ a~ i .i 
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Figure 62. Polarization versus scattering angle for dielectrics of Figure 61a. 
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Figure 64. Radiance versus scattering angLe for ice with distributions 

of particles for which a = 5, 15, 20, 25, and 30. 
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Thh radiance does increase slightly in both cases as a increases (Figs. 

max - 

64 and 67) . Due to the distribution, very few particles exist with a > 15. 

Their contribution to the total scattered radiance or polarization is small. 

The effect of cutting off the small particles for ice is seen in Figures 68 and 
69. The effect is especially noticeable in the polarization diagr&ms. If 
Figure 65 is compared to Figure 69, it may be seen that the particles with 

o; < 2 produce a strong positive polarization in the range 40 < 0 < 100 deg. 

This is caused by Rayleigh scattering from these particles. The dip at 
approximately 170 deg is due mainly to’ particles with size a s 5. 

Many cases occur for which a cloud of particles is not homogeneous. 
Results of radiance and polarization calculations for a half-and-half mixture 
of ice and quartz are given in Figures 70 and 71. As before, the backscatter 1 
peak of dielectric particles tends to indicate that these particles aye necessary 
for the zodiacal light. The positive polarizations of water in Figures 62 and 

65 indicate that dielectrics as well as metals could be candidates for the 
interplanetary material. 

. The normalized scattered intensities and polarizations obtained from 
Categories I through III distributions are shown in 72. A size range of 0. 1 
to 20 was used so that in the cases of Categories I and in both sides of the peak 
would be covered. Figures 72a through 72d show the radiance and polarization 
from the modified gamma distribution (specifically distribution 1, Table 4) ; 

n (a ) = (1.49 x 10 5 ) a 6 exp (-35.18 a °* 1705 ) 

This distribution was tailored to peak at a = 1.0 and to have a slope of -4 at 
oi = 20. The large particle contributions are evident in the large diffraction 
peak and also in the degree of polarization for ice. 

Figures 72e through 72h show the radiance and degree of polarization 
of ice and iron as obtained from an exponential size distribution, 

n(a) = (1.0 x 10“ 14 ) e 1 ^ 01 

The relative number of particles with size parameter a near 10 is larger for 
this distribution than for the particular Categories I and in distributions con- 
sidered. Curves for the degree of polarization have structure which indicates 
the effects of the large particles. The logarithmic distribution used is 
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Figure 66. Polarization versus scattering angle for iron with distributions 

for which O' „ = 5, 15, and 30. 
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SCATTERING ANGLE (DEG) 

Figure 67. Radiance versus scattering angle for iron in distributions 

for which n(o; ) ~ a " 4 and a = 5, 15, and 30. 
x ' max 
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Figure 68. Radiance of water in distributions for which 

a . = 2, 3, 4, and 10. 
mm 
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Figure 69. Polarization for water in distributions in which 

. =2,3, 4, and 10. 
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Figure 70. Polarization. of 50-50 mixture of water and quartz. 
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SCATTERING ANGLE (DEGREES) 

Figure 71. Radiance of 50-50 mixture of water and quartz. 
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Figure 72b. Normalized scattered intensity versus angle for ice in a 
modified gamma distribution (Category I). 
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Figure 72c. Degree of polarization of iron in a modified gamma 
- ■ ' distribution ( Category I) . 



Figure 72d. Degree of polarization for ice in a modified gamma 
distribution ( Category I) . 






Figare-72e.i Normalized scattered intensity versus angle for iron in an 
exponential size distribution (Category II). 



Figure 72f. Normalized scattered intensity versus angle for ice in an 
exponential size distribution (Category II). 
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Figure 72g, Degree of polarization of iron in an exponential 
distribution ( Category II) . 



Figure 72h. Degree of polarization for ice in an exponential 
distribution (Category n). 




n(a! ) = 0.7041 exp - 


/ in a \ 
\ 0.5nT2 ) 


2 


A modal size parameter of a = 1.0 and a distribution width of cr = 0. 5 were 
assumed. Kerker [ 24] points out that for a distribution such as this (shape 
factor m = 0) there is no structure in the radiance curves when u is as large 
as used here. For the distributions he considered, a a of 0.3 is the upper 
limit. For cr’s smaller than this, maximum and minimum peaks did appear. 
Figures 72i through 721 shows the radiance and polarization as obtained from 
this distribution. As can be seen from the smooth radiance and polarization 
curves, the position of the peak values of degree of polarization and the low 
values o/ radiance, the smaller particles dominate this distribution. The dip 
at, 175 deg in the radiance curve for ice does not appear. This dip is apparently 
due to particles with o; > 5. In this distribution, the number of particles with 
o' > 5 is very small. 


J. Ecliptic Zodiacal Light Models 


Mie theory has been used in models of the zodiacal light to calculate 
the components of the scattering intensity functions. Although the particles 
are restricted to spheroids, the models are on much firmer footing than those 
based on empirical laws. In a series of papers, Giese et al. [16, 26, 32, 35] 
considered the photometric properties of different models by the use of Mie 
theory. This section closely follows these papers. 


Start with the basic equation for the radiance of the zodiacal light 
[ equation (33)1 of the heliocentric model restricting the region of interest 
to the ecliptic plane. 



R E® 8, N 
X X o 

. m+1 
sin e 



where 




min 


(172) 
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Figure 72i. Normalized scattered intensity versus angle for iron with a 
logarithmic (shape factor m = 0) size distribution (Category in) . 



Figure 72j. Normalized scattered intensity versus angle for ice in a 
logarithmic (shape factor m = 0) size distribution (Category HI). 
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Figure 72k. Degree of polarization of iron in a logarithmic 
(shape factor m = 0) distribution (Category in). 



Figure 721. Degree of polarization for ice in a logarithmic distribution 
( shape factor m = 0) ( Category III) . 
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and 


N = / 

n J 


max 


mm 




da 


From these equations, one sees that the model parameters are: 

n = number density parameter, 
o 

a = size distribution parameter. 

a = maximum particle size, 
max 

- a . = minimum particle size, 
mm 

p = size distribution exponent, 
m = spatial distribution exponent. 

This general model assumes a density distribution of the form 


, , -m -p 

T)(a, r) = n Q r a 

A specific model calculation requires the selection of the refractive indices, 
or equivalently the selection of the material composition, and the percentages 
of each material selected. A series of selections and iterations will probably 
be needed before a model satisfying observations is obtained. 

Before considering different compositions, a single material (ice) will 
be considered and the above parameters varied. N q is varied to obtain a 

match to the magnitude of the observed radiance. The parameters a . , 

min 

a , p , and m are varied to study their effects on the elongation depend- 
max 

ence of the radiance and polarization. 

For a max small, all particles will scatter light approximately accord- 
ing to the Rayleigh law. The calculated radiance from a distribution of small 

particles (0.01 < a < O.l) with m= 1, p= 4.0, and N = 5 x 10 -2 particles/ 

o 

cm 3 is shown in Figure 73. The number density 
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a -p 

max 

N = — f n /— ) d a , (173) 

o 27r j . o l a I ’ 

a ■ \ o / 

min - 

is calculated at a wavelength of X = 0. 5 and a is set equal to 1. 

. o o 

The variation of radiance with elongation roughly follows the observed 
dependence with the obvious discrepancy at e = 180 deg. The polarization is 
shown in Figure 74. Since the scattering intensity functions for electrons 
behave according to the Rayleigh theory, this distribution of small particles is 
equivalent to approximately 1.1 x io i * 3 electrons /cm 3 , a value which is too high 
[57, 58] . The polarization is large for e< 90 deg. It peaks at 60 deg with a 
value of 0.60 which is three times greater than the observed value. 

In Figures 75 and 76, a is varied while a . = 1. For a particle 

max min 

distribution of a P , the radiance and polarization do not have a strong de- 
pendence on particles with a greater than 10 (Part I of this section) . In 

Figure 75, the radiance curve for a = 20 is very close to that for a = 

& max J max 

10. The curve labeled ’’coarse” uses a 16 -point Gaussian integration proce- 
dure which has been used in most cases in this study to save computer time. 
Although the coarse integration gives some false features, it is felt that the 
comparisons made are still valid. 

As one would have anticipated from the single particle and distribution 

results, the variation of a . with a fixed affects the degree of polariza- 

min max 

tion strongly, as shown in Figures 77 and 78 [22] . The increase of the 
radiance (Fig. 77) with increasing a is due to holding N q a constant and 

to the greater scattering efficiency per particle for the large particles. 

Figures 79 and 80 show the results of varying the spatial dependence > 
parameter m. The radiance increases for e < 50 deg as m increases and 
decreases for e > 50 deg as m increases. Polarization increases with m. 

From the earlier study of 

i n 

— — f sin m (0) d© , ' 

. m+1 J 
sin e e 
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Figure 74. Degree of polarization for very small (a = 0. 01 - 0. l) 
ice particles (N = 0.05 particles/cm 3 , m= 1, and p = 4) . 
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Figure 76. Degree of polarization for a 
variation of a for Figure 75. 
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Figure 78. Degree of polarization for variation 
of a . for Figure 77. 



ELONGATION 


Figure 79. 
(N = 1. 5X 10 -11 


Radiance for a variation of m 
, P = 4, a . = 1, and a 





it was noted that 


I (e= 7i) = — 

m v ' m + 1 



which shows the inverse dependence of the radiance on m at e = 180 deg; The 
increase of radiance toward the sun is caused by the increase in particle 
number (i.e. : , increasing m) . The increase of polarization with increasing 
m results from having less smoothing effects along the line of sight, since 
for higher m the contribution to the polarization by the larger scattering angles 
decreases. For comparison purposes, the results of Giese [26] and Giese and 
Siedentopf [ 16, 84] for the variation of m are shown in Figure 81. Their 
calculations are based on particle sizes that go from a = 2 to a = 120; however, 
the variational dependence is the same as shown in Figures 79 and 80. 

Table 24 compares the exponential gradient for the variations of m and 
p for the range of 30 to 60 deg., k is directly related to an increase in m and 
varies weakly with an increase in p. The variation with p is shown in Figures 

82 and 83. Increasing p gives more weight to the smaller particle sizes in 
the distribution. These particles have lower scattering efficiencies and a 
Rayleigh-like polarization. The decrease in scattering efficiency causes the 
radiance curves to be lower forhigher p. Likewise, the higher p, the polari- 
zation peak increases and shifts towards 90 deg when Rayleigh-type particles 
are included (cf. Ref. 22) . Giese’s results for varying p„are shown in Figure 
84. Except for the changes in magnitude of the radiance and polarization, I 

the dependence on p is the same. An explanation for the lack of changes in 
the shape of Giese’s curves for quartz and ice as compared to. Figures 82 and 

83 is that his distributions, 1 in starting from a = 6 and .a = 12, eliminate the 

particles which contribute most to the change. The structure ,of the curves 
for iron does not change because the scattering cross section after a « 5 is 
more stable than that for dielectrics. , 

TABLE 24. EXPONENTIAL GRADIENT (30 TO 60 DEG) FOR ICE 



K 


K 

m = 0 

2.6 

P= 3 . 

3.3 

1 

3.3 

4 ' 

1 

3.3 

2 

p = 4 

3.9 

5 

m = 1 

3.2 
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For a given m and p, the shape and magnitude of the radiance and 
polarization curves depend on the material. The results of varying the ma- 
terial while keeping m = 1, and p = 4 fixed are presented below. For these 
calculations, the particle size range is 1 ^ a ^ 20. 

In Figure 85 the radiance curves for the metallic materials iron 
(m* *= 1.27 - 1.31. i) , carbon (m* = 2.45 - 1.42 i) , nickel (m* = 1.46 - 2. 68 i) , 
and zinc (m* = 1.41 - 4. 10 i) are compared. Polarization and radiance for 
these metals have similar dependences on elongation, although the numerical 
values differ. No increase in the antisolar direction is seen in the radiance'. 

The polarization from 130 to 180 deg is very similar for all metals. Position 
and value of the maximum polarization for each material are given below 


(Fig: 86): 



•; Material 

Maximum Value 
of Polarization 

Position of 
Maximum (deg) 

,v ' \ Iron. 

" 0.49 

46 

Carbon 

. 0.30 

48 

Jr . - 

Nickel. 

0.31 

33 

Zinc 

0.24 

33 


The maximum in the polarization occurs in the forward direction. 

Dielectric materials are compared in Figure 87. Except for dirty ice 
(m* = 1.33 - 0.5 i) which approaches a metallic material, the radiance shows 
the effect of a gegenschein. The scattering efficiency increases as the 
particles become hard (i.e. , as the refractive index increases). Polarization 
fluctuates with no general trend (Fig. 88). Ice is dominantly positive, while 
quartz and titanite are dominantly negative for the range of particles con- 
sidered. 

Table 25 compares the exponential gradient for various materials. 
Generally the dielectric materials have a higher value of k than do the metallic 
materials for the 30- to 60-deg range. No conclusions about the composition 
of the interplanetary matter can be drawn from this because of the limited 
number of materials studied and the limited range of particle sizes used. 
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Figure 87. Radiance for a variation of dielectric materials 

(N = 1. 5 x lO -10 , p = 4, m = 1, a . =1, and a = 20) . 
v o mm max ' 
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Figure 88. Degree of polarization for a variation of metallic materials 
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TABLE 25. EXPONENTIAL GRADIENT FOR VARIOUS MATERIALS 
BETWEEN 30 to 60 DEG FOR THE MODELS CONSIDERED 



K 


K 

Titanite 

2.7 

Iron 

2.6 

Quartz 

3.0 

Carbon ’ 

2.6 

Ice 

3.3 

Nickel 

2.4 

Dirty Ice 

3.4 

Zinc 

2.3 


Since no single component model has resulted in duplicating completely 
the observed zodiacal light quantities [ 40] , various combinations of materials 
should be investigated. In general, metals give the correct position of the 
polarization peak, while dielectrics give the gegenschein effect. To combine 
these effects, a model should include at least one metallic and one dielectric 
material. 


For a combination of two materials, the radiance expression is 


L -I^ 1 ) + L^ 1 ) + L A < 2 > + L B < 2 > 
X 2 |_ X X X X 

« (j) 


7 A 7 B 
= L + L 
X X 


(174) 


where L^ VJ/ are the jth component and the a th particle type. The polariza- 
tion is 


l £U) + L B(1) . l A( 2 ) . l B(2) 
\ \ \ \ . 


P = 


(175) 


— B A 

In general, P is not equal to P + P , but the following statements can be 
\ \ \ 

made as a result of the definition above: 


( 1 ) if 


—A — r — 1 

L = Lf then P = - ( P“ + P' 
X X X 2 V X X 


/A B \ 

( P x + p x-) • 


(176) 
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(177) 


( 2 ) if 


<= p f <= p r= p * ■ 


(3) the polarization tends toward the polarization of the material with the larger 
total radiance. 

A method to fit experimental data using a two- component model is as 
follows (cf. References 16 and 84). From computed values of 


A(l) B(l) A(2) B( 2 ) 

L \ ’ \ * L \ * L \ 


(which could be curved fitted) , the values of 


7 A -B 


\ * V- p * 


B 


A B a B 

can be calculated for some N and N . Variation of N and N changes 

00 00 

“A — g A B 

the levels of L and L , while P and P remain unchanged. But from the 
\ \ \ \ 

statements above, the average polarization P will follow the material with the 

• A B * — 

larger radiance values. Hence, by varying N and N , a best fit for L and 

00 A 

P^ can be obtained. For a metallic -dielectric combination, the dependence 

of both the radiance and polarization is changed by varying N ^ and 

(i. e. , the number densities for each component) . The correct polarization 
peak is obtained from the metallic component, while the dielectric component 
gives the gegenschein or backscattering peak. Figure 89 illustrates the 

results that can be obtained by varying and for a given model. 

Clearly, almost any result can be obtained. The procedure can be extended 
to models with several compositions of materials. 


As an example, a dielectric-metallic and a dielectric-dielectric model 
are shown in Figures 90 and 91 along with single-material models. In the 

"• ) a B 

combination examples N = N (compare with Fig. 89). 
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Figure 91. Degree of polarization for combinations of quartz, ice, and iron. 



Giese and Siedentopf [ 16, 84] investigate combinations of electrons 
and purely dielectric dust; electrons, iron, and dielectric dust; and iron and 
dielectric without electrons. Their results are shown in Figures 92, 93, and 
94. The high electron values of Figures 92 and 93 result from early observa- 
tions which gave high values for the polarization. (Compare with Figure 74 
of electron polarization. ) 

Figure 95, which summarizes their results, illustrates that several 
composition models can yield approximately the same results, giving no 
uniqueness to the models [ 16, 84] . However spectral dependence was not 
calculated. Models I and II have £oo large a value for the electron density. 
Model in contains small (Rayleigh) scatterers. Weinberg [17, 18] concluded 
from his observations and the results of Giese and Siedentopf that a major 
fraction of both the radiance and polarization of the observed zodiacal light 
can be accounted for by small dielectric particles without resorting to high 
electron densities. 

Other composition models have been considered by other authors. 

Kovar and Kovar [ 85] based their compositions on the Apollo 11 lunar sample 
analysis. Consequently, they include glass in their models in varying amounts. 
The glass is considered to be in three forms having a refractive index of 
m* = 1, 55, 1.65, and 1.75. Several models were reported which differed in 
material^ and percentages of the constituents. Two models from their prelim- 
inary report are shown in Figures 96 and 97. The parameters for the models 


are: 
Model 5 

Percent 

Refractive 

Index 

m 

P 

Iron 

60 

1.27-1.371 

1.5 

2.5 

Ice 

30 

1.33 

0.0 

2.0 

Glass 

7 

1.55 

2.0 

4.0 


2 

1.65 

2.0 

4.0 


1 

1.75 

2.0 

4.0 
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Figure 92. Giese and Siedentopf’s [16, 84] results for 
various electron-dielectric models. 
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Figure 94. Giese and Siedentopf’s results for an isotropic reflecting 
dielectric and a hypothetical metal having 
Ingham’ s polarization function [161 . 
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Figure 95. Giese’s [26 ] results for radiance and polarization for three 
models assuming different electron densities. 
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Figure 97. Polarization for Kovar and Kovar's Models 5 and 7. 
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Model 7 

Percent 

Refractive 

Index 

m 

P 

Iron 

50 

1.27 

1.0 

2. 5 

Ice 

20 

1.30 

0.0 

2.0 

Glass 

20 

' 1.55 

0. 5 

3.0 


7 

1.65 

0.5 

4.0 


3 

1.75 

' 0. 5 " 

3.0 

N q = 5 x 10” 14 particles/cm 3 




Model 7 has less ice and more glass with a different spatial and size dependence 
than Model 5. The radiance curves for both Models 5 arid 7 are in good agree- 
ment with the observed radiances of Weinberg [47] and Wolstencroft and 
Rose [ 52-1. However, the polarizatiort.of. these models is somewhat lower < 
than the observed values; for e > 160 deg, the deviation is significant. 

Earlier, Little et al. [86] also considered a model with a glass content. 
Using a model with 97 percent ice (m* = 1.33) and 3-percent pyroxene (m* = 
1.70) , they were able to obtain moderate agreement with the polarization 
observations of Ingham and Weinberg. But, for e > 160 deg, their curves 
exhibit the same behavior as the curves of Kovar and Kovar given above. 

This is because of the ice content with a particle size range of 1 s a ^26 
and p = 4. Little concluded that for the polarization of this model to agree 
with observation, m must be greater than 0.6. 

Giese and Dziembowski [32, 33], in connection with out-of-the-ecliptic 
analysis, give an' iron-quartz model and an ice-dirty ice model which agree 
approximately with observations. Their models (Table 26) are based on 
scattering by absorbing and dielectric particles large enough to produce non- 
Rayleigh scattering. The different models use the fact that the scattering 
functions for a distribution are dependent on the imaginary part of the refrac- 
tive index. For example in the distributions used by Giese and Dziembowski, 
the polarization for particles with m* = 1. 0 - 0. 01 i peaks at — 90 deg; for 
particles with m* =1.0-2.0iit peaks at « 40 deg. Both models used by 
them have a particle size range of 2 to 120 and a size distribution with p = 2.5. 
They noted [33] , models which contain large particles, whose sizes are 
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distributed according to the inverse power law, allow an increase of the 
particle sizes to increase the forward portion of the scattering intensity 
function without seriously altering the polarization of higher scattering angles. 
This differs from the case with very small particles (a, = l) where any con- 
siderable increase of particles with a high a value would result in deviations 
from the Rayleigh-like scattering behavior. 

To match observational data, using single material models, Aller et al. 

[ 22] had to ’’break 11 the power law distribution (Table 26) ; i, e. , for some range, 


a . 
min 



P = P L 


and for 


(178) 


a < a ^ a 
B max 


P = P 


M 


(179) 


Any number of breaks could (be included, but they employed two break points. 
Using dirty ice (m* = 1.33 - O.Ol) , m = 1 and 


1.25 — a — 25.5, 

they obtained an approximate agreement with observations in three cases; for 


a . = 0. Ol/i 
min 


P-_ = 4. 0 
M 


and 

P l = 2.0 
P L = 2.3 
P l =3.0 

However, in breaking the distribution, another parameter, the break 
point radius a^, is introduced into the calculations, as well as another ex- 
ponential parameter for the distribution. 


a B = °*40 

a B = °* 50 

a = 0. 60 
13 


The same situation is obtained with the belt models. The belt models 
increase the versatility of the model but also increase the number of parameters. 
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Powell et al. [28] matches the observations of Weinberg [46, 47, 48], 
Blackwell and Ingham [9, 10, 11], Divari and Asaad [64] ,-,-Robley [45], and 
Behr and Siedentopf [8l using such models.! They conclude that particle sizes 
are in the range of a = 0. 08 to 1. 5. The number density of particles in each 
size range is approximately constant between the orbits of Venus and Mars 
where dielectric particles (l.4=< m < 1.8) dominate. Powell et al. concluded 
that the density of particles in the region beyond Mars is not sufficient to 
affect the character of the zodiacal light. 

In these models, the size distribution is arbitrarily truncated at some 

lower limit a . . Although this is within the order of the current models, it 
min 

is not physically realistic, especially in certain cases (appendix). 

In summary, many models can closely match the observational data 
taken in the ecliptic plane. Although good physical arguments may exist to 
justify preference of some models over others (e.g. , multimaterial versus 
singl e-material models) , these arguments have no firm basis. Table 26 is 
a partial listing of zodiacal light models which are in close agreement with 
observational data. .... • 

K. Spectral Dependence and Null Point 

An additional step in studying the behavior of the radiance and polariza- 
tion from models is to study the variation of these quantities with wavelength. 
Color differences between models is a function of the refractive indices of the 
material used. Thus, observations of radiance and polarization in several 
wavelengths would help establish possible compositions of the zodiacal cloud. 


Color magnitudes are calculated by 

M a = - 2. 5 log 10 (L^) . . .1 (180) 


For two wavelengths, the color excess, with reference to the spectral content 
of the sun, is calculated from 


A (M A - M b ) = (M a - M b ) + 2. 5 log 10 


^A 

*\B 


(181) 
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L . is the radiant energy flux, at wavelength X., from the sun. Color excesses 

Al 1 ... t 

-for Rayleigh scattering are easily calculated. Since the scattered energy flux 
is » the color magnitudes are given by 

' ■■ 

M a = - 2. 5 log 10 -jr -2. 5log 10 R . 

All terms common to each calculation are included in R. The color excess 
is thus 

A(M a -M b ) = - 2. 5 log 10 + 2. 5 logi 0 jr- 

■ B 


+ 2.5 log 10 


^XA 


I =10 log 10 \ A - 10 log 10 \ B . 


‘XB 


Using explicit subscripts and the wavelengths \_. = 380 Nm, X = 500 Nm, 

V (jr 

and X = 650 Nm, the color excesses 

XV - . * 


A (M y - M Q ) = -1.192 , . 

and .... ■ ... .*-■ ' • 

■' ■ , A (M - M ) = -1.139 

XV * , r • . t 

are obtained. Thus, the color enhancement in both cases is bluish. These 
values represent an upper limit for small particle models. Aller et al. [22] 
using these same wavelengths, calculated the color excesses from Model FI 
of Table 26. Selected values of their results are listed in Table 27. From 
their Mie theory models they concluded that the zodiacal light is about 0.2 
magnitudes bluer than the sun in both V-G and G-R color indices. Powell, 
et al. [28] , using the belt model, also found some color enhancement in the 
blue in agreement with observations of Tanabe and Huruhata [ 87 ] . : Blue 
enhancement indicates submicron-size particles because of Rayleigh scattering 
Gindilis and Karyagina [88] find a wavelength dependence of the zodiacal light 
that goes as X -0 * 1 , while for the gegenschein the dependence is X 1 * 28 ; 
Peterson [43, 893, on the other hand, finds a x 0,065 dependence, a reddening 
effect. (Peterson used a G1 star for his calculations.) 
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TABLE 27. COLOR EXCESSES FROM ZODIACAL LIGHT MODEL Fi [22] 


e 

A (M v -M q ) 

a(m g -m r ) 

Degree of 
Polarization 

Approximate 

Radiance a 

20 

-0.132 

-0.252 

0.063 

50 000 

40 

-0.152 

-0.061 

0. 129 

6 310 

60 

-0.148 

-0. 170 

0.218 

1 590 

80 

-0. 197 

-0.157 

0.236 

631 

100 

-0.203 

-0.181 

0.213 

399 

120 

-0.187 

-0. 228 

0. 102 

317 

140 

-0.273 

-0.211 

: -0.061 

282 

160 

, -0.054 

-0.239 

-0.068 

317 

i 


, a. Read from Figure 11 of Reference 22. 

The dependence of the neutral point position of polarization on wave- 
length is not well defined because of a lack of observational data. The de- ’ 
pendence found by Weinberg [40] has yet to be predicted by any published model. 
For ice (n ~ a ~ 4 ) , the null point variation is not that which is observed. Fig- 
ure 31 shows the null point variation for our basic model (1 ^ a ^ 20, n- a" 4 , 
r" 1 , ice). The null points were determined by extrapolating from the values 
at e = 150 and 160 deg. The tendency is for the null point to move' to higher 
elongations with increasing wavelength rather than to decrease as observed. 
Greenberg [90, 91] finds this is the behavior of the null point for ice. How- 
ever, for particles with m= 1.65 whose sizes are distributed according to 

n(a) ~ e" 5 ^ 43 ^ 3 , he finds, over a limited wavelength region, that the null 
point position moves to lower elongations with increasing wavelength as 
observed. 

• \ . 

Greenberg has observed that mixtures of metallic and dielectric materials 
(e.g. , Giese’-s models) should be examined. A significant difference between 
them would be expected because the backscattering from metals differs from 
that of dielectrics. 
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L. Out-of-the-Ecliptic Model 

. . . •- • ■ • - - - 

The out-of-the-ecliptic models discussed in Chapter in were: 

1. Constant thickness model [27] . 

2. Linear thickness model [27] . 

3. Ellipsoidal symmetry model [32, 33]. 

4. Exponential model [22], 

5. Meteor model [35, 36]. 

i ' • 

This section presents the results of radiance and polarization calcula- 
tions using the models given above. Only those values of the parameters 
which made the model fit the observations closely will be considered. How- 
ever, observational results out of the ecliptic are scarce. 

j 

Gillett [27] used an empirical formula for the total differential cross 
section do^/dSl (Fig. 39) in the constant and linear thickness models. The 

specific cross sections used and the values obtained for K. and (2R/W) for the 

constant thickness model and K. and (u> ) . for the linear thickness model are 

given in Table 28. W is the width of the interplanetary cloud, u> is a defining 

cone angle, and K. is a weighing factor (Fig. 17). The out-of-the-ecliptic 

polarization obtained by the models was high compared to that observed. Gillett 
obtained results for the constant thickness model with m = 1/2 or 1 and for . 
the linear thickness model with m = 1/2 which were fairly consistent with 
observations in the ecliptic. The variation of L (45, i)/L (45, 0) for various 
values of (2R/W) is shown in Figure 98. The best fit is obtained by having 
da^/dfi reproduce L (e, 0) for 30 to 100 deg with L (90, 90) = 0, such that 

L(90) > L(180) > L(110) with the minimum value L(e) in the range 125 < e < 
155. 

To under stand the behavior of Gillett-type models, several simple 
models and several models using the Mie theory were considered by the 
authors., Therresults are shown in Figures 99 through 109. In Figures 99 

/J 

and 100, the radiance contours for a hemisphere are shown for a simple 
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GILLETT’S RESULTS [27 ] 
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L 71 ;(45°, i)/L yi <45°, 0) 


m = 1 


(1) 2R/W = 2.5 

(2) 2R/W = 5.0 

(3) 2R/W = 7.5 

(4) 2R/W = 10.0 

(5) .5v(2R/W = 2.5) + .5 -,(2R/W = 7.5) 

AND OBS 




ECLIPTIC INCLINATION 


Figure 98. Gillett's [ 27 ] results for a constant thickness model. 
(The variation out of the ecliptic for the radiance is shown 
for various thickness models. ) 





M = 0 CONSTANT THICKNESS 



Figure 99. Out-of-the-ecliptic results using the constant thickness 
model and the basic heliocentric model (R/W = 2.5) 


M » 0 ' 

LINEAR THICKNESS 


90 ° 



Figure 100. Out-of-the-ecliptic results using the linear thickness 
model and the basic heliocentric model. 


heliocentric model with constant thickness and with linear thickness. The 
ecliptic pole is at the center of the figure and the ecliptic plane is on the 
circumference. The difference in the contours is that the linear thickness 
model is only piece-wise smooth and has a discontinuity in the rate of change 
of the radiance when the line of sight reaches the defining cone. The constant 
thickness model has the minimum radiance shifted from the ecliptic pole, whereas 
it is at the pole for the linear thickness model. Changes in the radiance and 
polarization of the constant thickness model with a change of R/W (=ROW) 
are shown in Figures 101 through 104. Ice (m = 1.33) was used as the particle 
material. As R/W is decreased (i.e., W is increased), the* radiance values 
off the ecliptic increase although the values on the ecliptic remain the same 
and the minimum value remains in the same region. The value of R/W = 1/10 
is the best fit to radiance-observed data. The polarization value increases 
for e < 90 deg (Fig. 104) and decreases for e > 90 deg. For this figure, the 
point of observation moves such that e = 90 deg is the ecliptic pole. -Figures 
105 through 107 show the result of changing the materials. The contours are, 
in general, established by the distribution scattering functions, as is the 
effect of changing p (Figs. 108 and 109). Refer to Figure 102 for a compari- 
son. 

The ellipsoidal symmetry model [32, 33] adopted a spatial distribution 
out-of-the-ecliptic mode that is continuous, ; 



Two of their models (Table 29 and Fig. 110) approximate the out-of-the- 
ecliptic results of Smith, et al. [19], In Figure 110, St is the ecliptic latitude 
and 6 is the distance from thie sun in ecliptic longitude. 


TABLE 29. GIESE-DZIEMBOWSKI MODELS [33] 


! • 

Model 

Refractive 

Index 

Alpha 

Range 

N /particles'^ 
o \ cm 3 / 

m 

{ 

r 

, I 

1.5 

6 to 120 

5.3 x iO" 15 

3.0 

10 

Model 

1. 5 

6 to 120 

4. 2 x 10 ,1 

0.1 

10 


. 1.27 - 1.37 i 

2 to 120 

2.1xl<f i3 

3.0 

10 

n 

1.33 

12 to 120 

8. 7 x 10-J® 

0.1 

300 

Model 

1.33 - 0.1 i 

2 to 120 

7.5xl0* 13 

0.1 . 

300 
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Figure 110. Giese and Dziembowski’s [33] results for two ellipsoidal 
symmetry models compared with Smith et al. [ 19] observations. 

Model I is metallic and Model n is dielectric. Both models have large 
particles with slowly decreasing (p = 2. 5) size distribution. By comparing 
large particle models with Rayleigh models (i.e. , smaller size particles) 
which tend to have a minimum further from the sun than the ecliptic pole, 
Giese and Dziembowski concluded that the larger particles better fit the ob- 
served data. 

Van Rhijn used a density distribution of the form 


n(x, y, z) = N q (1 - a (x 2 + y 2 ) - bz 2 ) , 


( 182 ) 
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where x, y, z are heliocentric coordinates for out-of-the-ecliptic calculations 
[ 2] . The integral was terminated when 

n(x, y, z) = 0, 

and a phase function was assumed; thus, 

4n h(0) = 0.655 - 0.711 sin 6 - 0.345 cos 0 + 0.406 sin 3 0 (183) 

for the angular distribution of the . scattered light. The values of the parameters 
a and b were determined from observed radiance to be a = 0. 176 and b = 1.406. 
This gives a ratio of the major-to-minor axis of 2. 8. For geocentric longi- 
tudes (ecliptic longitudes) of 6=0, 20, and 60 deg, the values of Van Rhijn 
are compared with the out-of-the-ecliptic data of Smith et al. [ 19] in Figure 111. 
The data points (6=0 and e = 80 deg) were used to normalize the data. The 
agreement is good. Using a particle with an albedo of 0.3 and a 50-cm radius 
(rocks) , Van Rhijn estimated the density of particles at 1 AU to be 10 -18 
g/cm 3 . 


Fesenkov [ 34] also used the ellipsoidal model and the scattering func- 
tions for aerosols to determine the polarization at the ecliptic pole. The 
following values resulted: 

. Eccentricity Polarization 

e = 0 disc model p = 34. 0 percent 

e = 0. 1 p = 32. 8 percent 

e = 1 (spherical model) p = 26.4 percent 

Aller et al. [ 22] used the exponential model to calculate out-of-the- 
ecliptic contours. However, their polarizations are generally too high com- 
pared with observations; cf. , Table 13, p. 107. The radiance minimum 
occurred near e = 130 deg, with i = 90 deg. However, the exponential decay 

_5z/ R 

away from the ecliptic varies as e , which gives an exponential half- 

width of the cloud of on the order of 10 1 AU. 

Rouy and Aller [ 92] make the observation that the distribution of the 
radiance perpendicular to the ecliptic could be approximated by a Gaussian 
distribution. For scans orthogonal to the ecliptic plane at a specific elonga- 
tion e, the radiance could be expressed by 
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ELONGATION 


Figure 111. Van Rhijn’s [3] ellipsoid results for out -of-the -ecliptic 
(solid lines) compared with Smith et al. [19] data (• ~ 6 = 0 deg, 
O ~ 6 = 20 deg, □ ~ <5 = 70 deg) normalized at 
6=0 and e = 80 deg. 





L (oj, e) = L zq (e) e ^ + Lg (c) ' 

where ■ 

L (e) = the peak radiance of the zodiacal light, 
zo 

L (e) = the average radiance of the background. 

13 

t 

oo = the angular measurement from the ecliptic. 

I x 

By analysis of the derivatives j various zeroes (Fig. 112) arise which one can 
use for determining the parameters in the equation. Examining Weinberg’ s 
data for a particular scan for total radiance, R’ouy and Aller determine that 
K = 1. 182 (Fig. 112) . Their calculation should be considered a preliminary 
one since the accuracy of the approximation has yet to be determined. , 

Singer and Bandermann [ 21] used the generalized model with the size 
distribution given by 


00 


n(r,9) = - J f f 


7t/2 


r/2 | r/a-1 1 tt/2-0 


a 2 f(a, e, i) da de di 

[ e 2 - (r/a-1) 2 ] (sin 2 i - cos 2 0 )V 2 

■ ' ' (184) 


where 


-bi 


f(i) = sin i e 

They assumed a scattering intensity function i(0 ) of the form 


V e > ■ 2 


(185) 


/ 1+T7 2 \ rj© 1 y < S > x _3 / v 

(— ^r) e • < 186 > 


where the first term is the diffraction component, and the second is the re- 
flection component, y is the albedo of the particle, -q is the phase coefficient of 
of reflection, and <S> and <S 2 > are the mean and mean-square of the particle 
size dimension. The values ofrj = 0.02, T= y < S >/< S 2 > = 0.14, b= 3, 
y = 0.1, and m = 1. 5 fit the observations of Smith et al. [19] best for the 
functions: 
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and 


L(e)/L(e = 30 deg), i= 0 deg 
L(i)/L(i= 0 deg), e = 90 deg 


L(e, i = 0 deg)/L(e, i = 90 deg) . 

T = 0. 14 indicates a flat size spectrum of large particles and b = 3 indicates 
an average orbital inclination of i= 30 deg. The radial dependence is r ‘ . 
Polarization values were not calculated because of the form of the scattering 
intensity function, 

Divari [ 93] has taken the meteor model with a scattering function of 
the form 


f (0) = 0.0549 [ 1 + 11.1 (e - 0.009)] (187) 

and calculated relative radiance. The isophotes obtained are close to those 
observed. The polarization was also matched by using an empirical function 
taken from atmosphere observations. The zodiacal light polarization values 
obtained are high, even in the ecliptic (e.g. , p(e = 80) = 40 percent); hence, 
due to the values of the parameter used as input by Divari, the results are 
inconclusive. 

It can now be seen that there are several out-of-the-ecliptic models 
(some using Mie theory, some not) which achieve reasonable agreement with 
observations. The agreement is not yet close enough to be realty satisfactory, 
but further improvement will probably await the results of measurement made 
outside the earth’s atmosphere. 

M. Nonspherical Particles 

The Mie-theory models deal only with spherical particles and can show 
peculiar effects, e.g. ; haze bow [ 32] . Analysis using this theory may give an 
erroneous interpretation of experimental data, especially in the anti-solar 
direction where the geometry of the particle plays a more important role [ 90, 
91] . However, Little et al. [86], Napper and Ottewill [ 94] , and Donn and 
Powell [ 95] have shown that the radiance, in the 50 < 9 <150 deg range, from 
a distribution of nonspherical particles of random orientation can be represen- 
ted approximately by an appropriate distribution of spherical particles. For 
scattering angles 0 < 60 deg, the protrusion of irregular structure increases 
the scattering function by diffraction. Nonspherical particles may prefer a 


248 



given orientation, resulting in further differences from Mie theory. Though 
Powell et al. [28] concluded that the scattering intensity functions at a given 
wavelength for a distribution of cubes (large- volume shape) could be repro- 
duced by spheres whose distribution is approximately that of the cubes. They 
also conclude that a size distribution for small- volume shapes may differ 
markedly from that of equivalent spheres. The difference between small- 
volume scatterers and spheres can be seen geometrically by considering, the 
ratio of the projected area of rods and disc to that of a sphere. If the rods 
have a length-to-radius ratio of 10, and the disc has a radius thickness ratio 
of 10, then, for an equivalent volume, the projected area ratio to a sphere is 
1.7 and 5.5, respectively. Hence, small-volume particles, such as flat 
particles and very elongated particles, are more efficient scatterers than 
spheres for equivalent volumes. 

Only a few scattering problems for finite bodies have been solved 
{e. g. , sphere, prolate spheroid, oblate spheroid, disc and rods [96]} although 
various perturbation procedures have been applied to other shapes. Recent 
investigations on finite cylinders might lend themselves to zodiacal light 
studies [24, 97, 98, 117, 118, 101] . Cylindrical scattering differs from 
spherical scattering in that there is a loss in the degree of symmetry which 
introduces anisotropic scattering. A shape factor can give a general charac- 
terization of the anisotropic scattering. The anisotropic processes should 
lead to new results in zodiacal light model parameters, although no model 
exists employing cylinders to date. 

Figure 113 and 114 give the scattering intensity functions and polariza- 
tion for infinite cylinders whose axes are perpendicular to the incident direction 
[82, 98] . The polarization for the cylinders is not zero in the forward or 
backward direction, a marked difference from spheres. For larger a , the 
diffraction peak (Fig. 113) is observed as an increase in the radiance in the 
forward direction for all angles of incidence. Kerker et al. [24] give single- 
size scattering functions for oblique incidence. The functions are strongly 
dependent on the orientation of the cylinder axis with respect to the incident 
direction. The scattered intensity functions contain a perpendicular and a 
parallel component for both perpendicular and parallel polarization of the 
incident beam. 

For a distribution of dielectric cylinders, (m* = 1.33) , Hanner [97] 
observed that with natural light at normal incidence, the total scattering 
intensity is analogous to dielectric spheres, having a forward diffraction peak 
and an increase in the back-scattered radiance. However, the polarization 
differed significantly from that for spheres (Fig. 115) . 
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Figure 115. Hanner’s [97 3 results for cylinders with a distribution 
of the form n = n exp l -5(a/ao) 3 ] , a = tilt angle. 











For a distribution of particle sizes, as the tilt angle increases, the 
total intensity decreases. The tilt angle is the angle in the incident plane 
between the incident direction and the perpendicular to the cylinder axis. 

Also, as the tilt angle increases, the polarization in the forward direction 
does not change much while the polarization in the backward direction oscillates 
in polarity. Compare Figures 115a and 115b [97] . Making use of an exten- 
sion given by Lind [99] to the infinite cylinder treatment, Hanner gives results 
(Fig. 115c) for randomly oriented long, finite dielectric cylinders. The shape 
of the curve is similar to that for a distribution of dielectric spheres [m* = 
1.33, n(o;) = n 0 aT 4 , l<a< 20] (cf. Fig. 62). The differences are that 
the polarization for cylinders is nonzero in the forward direction, less nega- 
tive in the backward direction, and its peak is at 110 deg instead of 120 deg. 

Thus, it is seen that toe radiance and polarization from models depend 
greatly on the shape of the particles. Hence, the Mie theory approach can be 
considered a first-order approximation. To obtain better models, the effects 
of irregular shapes on toe radiance, polarization, and color of zodiacal light 
require further investigation, both experimentally and theoretically, and more 
reliable and consistent observations are needed. 

N. Areas for Further Research 

Presently, a unique model does not exist which predicts the observed 
photometric properties of toe zodiacal light. However, comparing the various 
models which are in approximate agreement with the observations (cf. Table 
26) , one finds that the following properties are representative of the zodiacal 
light dust cloud: 

1. Composition — Largely dielectric particles with a refractive 
index in the range 1.3 to 1.8 with smaller percentages of metallic particles. 

2. Shape structure — An equivalent size distribution of spheres in 
the range from 0. 1 to 1.0' micron or 1. 0 to 10 microns. Assuming the smaller 
size distribution, the exponential power law exponent is on the order of 3-4, 
while the larger size distribution has a smaller value of, say, 2-3. 

3. Spatial structure — The spatial structure of the cloud can be 
considered to be heliocentric with a decrease in the particle density according 

-m 

to r wher;e 0 < m < 3. 0. The out-of-toe-ecliptic spatial density models 
indicate a characteristic half-thickness of 10" 1 AU for the dust cloud. The 
number density is on toe order of 10“ 13 - 10 -15 particles/cc. 
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In conclusion, it can be said that an enormous amount of work has been 
done in relating light scattering models to the photometric observations of 
the interplanetary dust cloud; however, a conclusive knowledge of the physical 
properties of the dust cloud still remains elusive. 

The following areas need refinement for model calculations and 
analysis of the zodiacal light by photometry: 

1. The size and spatial distribution needs to be based on experi- 
mental data. This is particularly true in the case of particle-size distributions 
where a lower cut-off limit is assumed. A sharp cut-off seems unlikely. 

2. A set of polynomials needs to be created to form do/dfi . 

3. The effects of irregularly shaped particles compared to spheres 
should be better defined to categorize differences. 

______ 4. The methods of inversion of the photometric observations of ' - 

da/dS2 to find particle sizes and composition need more definition and study. 

5. The complete effects of the physical forces on the dynamics of 
the interplanetary particles need to be determined. 

6. If possible, a Mie theory model which satisfies all the gross 
properties of the observational data should be found. 
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APPENDIX. DYNAMICAL FORCES 

This appendix considers the physical' forces which act oh the interplane- 
tary particles, thereby possibly determining the spatial and size distribution 
and the temporal variation. 

1. Radiation Pressure [24] 

According to Maxwell’ s theory of electromagnetic radiation, light 
carries momentum as well as energy. The direction of the momentum is in 
the same direction as the energy flow and is given by 


M = E/c 


where E is the energy and c is the velocity of light. Electromagnetic momen- 
tum produces a radiation pressure on the interplanetary particles. The 
momentum of the incident field is scattered and absorbed by the scattering 
center. The total energy per second removed from the incident beam of unit 
radiance by absorption is equal to the cross section for extinction, C . How- 

ever, of the momentum carried away by the scattered radiation, that part which 
is associated with the forward component is restored to the incident beam. 

The time- averaged radiation force is [ 24] - ' - . c. 



r 


Z^7ra 2 

c 


( C - cos 6 C ) 
\ ext sea/ 


where 2^ is the spectral radiant flux density at the surface of the particle, 
Csca is the cross section for scattering, and 

gj" / + 1 ( ii + i 2 ) cos 9 d ( cos 9 ) 

cos Q = 

7 r z 1 -f j . . . . 

f - J-i + y d ( cos 6 > 



+ i ? ) cos 6 d(cos 9) / C 
v sea 


The asymmetry factor, cos 6 , is the mean of cos 6 with the scattering inten- 
sity functions as the weighting function. The efficiency factor for radiation 
pressure is defined as 




cos 9 Q 


A 

sea 


ext 


cos 9 


C 

sea 

7ra 2 
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The asymmetry function can be expressed as 


CO 


cos 6 Q 


sea 


^2 

n=l 


n ( D + .. 2 ) r 6 ( a * a + b * b ) 
n + 1 \n n + 1 n n + 1/ 


( 2n + l) „ 

+ , 4- Rea *b 

n ( n + 1) n n 


Hence, for a specific wavelength the pressure exerted on a particle of a cross- 
sectional area 7ra 2 is 


P X = —L = q X 

7ra 2 c pr 


2. Gravitational Force Versus Radiation Force 

The force of gravity on a particle of radius a, mass density p , and 
distance r from the sun is 

GM , 

_ © 4 a 

F = ~5 — — 


g r z 3 


where G is the universal gravitational constant ( 6. 670 x 10 11 nt m 2 /kg 2 ) 
and M is the mass o: 

is) 

and a in meters [ 100] 


and M is the mass of the sun ( 1. 989 x 10 3 ° kg) . For r in AU, p in kg/m 3 , 

is) 


2.5 X 10 2 , r „ 

F = 5 P a [nt] • 

g r 2 

To compare this to radiation pressure, one sees that the radiation force is 

F = ~2 Q HO R 2 f 2 A dX = — ^ f Q X (m, a) 2 X dX , 

r cr 2 pr ■ J c r 2 J ' 


256 



where 2 is evaluated at 1 AU. Q (a") is a function which is on the order of 

> pr' 

unity for a greater than 0.6 x 10 7 m and which declines approximately as a 4 
for smaller a. For r in AU and a in meters, 


F 

r 


-5 a 2 — 

1.4 x 10 y ~ Q [ntj . 
r 3 pr 


The quantity of 1 + p, defined as the ratio^of the force caused by radiation 

pressure F to the force of gravitational attraction F , is independent of the 
r g 

distance r of the particle from the sun and is given by Gindilis et al. [105] 

F ^-2 ( \ \ 

1 + p = — = — Q (m, a) 2 d\ , 

F pa J pr 
g *i 

where p is the density of the dust particle, A = 3 R 2 /4cG . They use 
$ (m, a) = p (l + p) 


to represent the ratio of radiation force to the gravitational force per unit 
density for a given particle with radius a and refractive index m*. This 
quantity is important in evaluating the relative magnitude of the two forces. 
Table A-l gives the values of Q (m, a) obtained by them; Table A -2 gives 

the values $ (m, a) ; Figure A-l gives the pertinent results for several 
substances. In the case of quartz and water, the radiation pressure is less 
than the force caused by gravitational attraction for all particle sizes. 


An analogous situation occurs for most dielectrics (m* s 1.6 and p 
~ 2-3 g/cm 3 ). When m is larger than two, particles of certain sizes will 
be pushed out of the solar system. For titanite (m* = 2, p = 3.5 g/cm 3 ) , 
particles with 0.1 s as 0.3* 10 -6 m will be ejected. If Q ^(a) = 1, then the 


_ n 

radiation pressure dominates solar gravity for a s 5 x 10 m. The rapid 
decline of Q (a) below 1 x io -7 m limits the region of radiation pressure. 
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TABLE A -1. EFFICIENCY FACTOR Q (m,"a) [105] 
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TABLE A-2. VALUES OF $ (m. a) OR 1 + m WITH p = lg • cm' J [105] 



























Figure A-l. The magnitude of 1 + p as a function of the particle radius: 

1 = water, 2 = quartz (m* =1.5, p = 2. 5) , 3 = titanite (m* = 2.0, p = 3.5) , 
4 = iron (p = 7. 8) , 5 = graphite (p = 2. 0) [105]. 
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3. Lorentz Force 


The processes which are important in determining the electric charge 
q on an interplanetary dust particle are shown in Figure A -2. 

Because of magnetic fields carried in the solar wind, the particle can 
experience the Lorentz force which is exerted by the magnetic field B on the 
particle moving with a velocity v and is given by 


F = - qv x B 
L 

For I v | = 400 km/s and I B I = 3 x 10 5 gauss at 1 AU, then with a in 
meters [99] 


F « 10 a[ nt] . 

L ' 

The Lorentz force F is caused by the interaction with the interplanetary 

magnetic field alone, whereas the convective force F con " qw x B depends on 

the fact that the solar wind carries the magnetic field along at a speed w, which 
is large compared with v. 


4. Corpuscular Drag 

The drag forces caused by collisions with solar wind particles can be 
calculated by: 

F , = rate of momentum transfer 
cd 

= proton flux x momentum x area of particle. 

— 19 

An estimate of the total drag force is 6. 3 x 10 a 2 [ntl . The component of 
drag parallel to the velocity is 4 x 10 _2 ° a 2 [nt] [77] . 
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REFLECTED ELECTRONS 



Figure A-2. Processes which determine the electric charge 
on an interplanetary dust particle [102] . : 


5. Coulomb Drag 


An equilibrium charge q on a dust particle can be approximately calcu- 
lated if the charging process is the photoelectric effect and the discharging pro- 
cess is the collection of electrons from the solar wind. An estimate of the 
equilibrium potential is ^ 2 volts [ 103] . If the charged particle sees an elec- 
tric Field, then the particle will experience a Coulomb drag force, F = qE. 

-19 ^ 

An estimate of this force is 10 a 2 [ nt] and the drag component parallel to the 

particle velocity is 6. 6 x 10~ 21 a 2 [nt] [27]. 
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6. Mass Loss Effects 


The solar wind will cause a radial sputtering erosion. Solar radiation 
will cause evaporation or sublimation. Both of these effects contribute to mass 
loss of the interplanetary particles and are more important closer to the sun. 
Mass loss effects .will affect the dynamics through the equation of motion. 


7. Physical Process 


Sputtering by protons and alpha particles is most effective at kilovolt 
energies. The threshold energy for sputtering of metals by ions has been 
found to be given by the requirement that the product of the momentum trans- 
ferred at .threshold from an ion to the target atom and of the sound velocity in 
the target be proportional to the heat of sublimation of the target [ 36] . The 
erosion rate of a dust particle is proportional to the flux of ions in the solar 
wind. The mass rate loss is 


dm 
dt : 



> 


where da/dt is the shrinkage rate: 


a = radius 


da 

dt 


c m X F. Y. 
t u li 
i 


and m is the atomic mass of the target, F. represents the various fluxes, and 

V , ‘ • 1 

Y the corresponding yields ( atomic/ion) . 

To estimate the temperature, the solar irradiance blackbody equilib- 
rium temperature is calculated from the equilibrium equation; energy received 
is energy radiated: 

A. = ecr T 4 A 
r^ 7r t 

( Stefan- Boltzmann Law) where 

a = absorptivity 

S = solar constant [1390 watts (AU) 2 /m 2 ] 
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r = solar distance ( 1 AU = 1.496 x 10 11 m) 

A = projected area (solar-oriented) = Tra 2 (sphere) 
e = emissivity 

a - Stefan- Boltzmann constant, 

T = equilibrium temperature, 

A = total area = 47rr 2 ( sphere) 
t 

and heat loss from vaporization is neglected. 

Assuming that ( l) the sphere is rotating or the sphere is small enough 
that equilibrium temperature has meaning and (2) a = e; then for interplane- 
tary spherical blackbody particles the equation for the equilibrium temperature 
for various solar distance is 


T 2 r = 0.785 x 10 5 l°K 2 AU] 


8. Poynting-Robertson Effect 


Because light has momentum, its reradiation or general reflection 
from a particle in orbit about the sun produces a retarding force on the motion 
of the particle. While the process of absorption and re-emission produces no 
net force when one chooses to work with a stationary frame referred to the 
particle, it is found when the solar reference frame is used that a resisting 
force is introduced on the particle which is proportional to its velocity [ 104] . 
Assuming (1) spherical particles of radius a and uniform density, (2) gravi- 
tational forces are dominant over the radiation pressure, (3) diffraction is 
not considered, and (4) the particles absorb all incident radiation over a 
cross section 7ra 2 and re-emit the radiation isotropically at the same rate, 
then the equations of motion are 



J 1_ 

r 2 


2a r 
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where 


M = GM 0 - ac 

O' =3 E' 0 / 16 7TC 2 ap 

E’ 0 = total energy emitted by the sun per second (3.79 x 10 33 erg/s) . 

. • '■ \ 

The resultant equations for the semimajor axis b and the eccentricity e are 

' i 

db _ a ( 2 + 3 e 2 ) 1 


de _ 5a e 

dt 2b 2 (l-e 2 ) /2 

For circular orbits, the total time of fall for a particle of radius a, density p 
at an initial distance r in AU is [ 119] 

• i 

b 2 

t = — — = 7. 0 x 10 6 a pr 2 [years] 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, March 30, 1972 
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